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__ The Joint | Committee has approved the | publication of this Commentary on 
- Plastic Design in Steel with a strong ‘recommendation that it be read carefully — 
_ by all structural engineers. It is based ona series of preliminary reports 
prepared for the Welding R Research Council and the American Society of Civil - 
‘Engineers by a research group at Fritz Engineering Laboratory, Department — a 


ora _ Although much of the experimental and theoretical work reported was per- — 
- formed at Lehigh University, the WRC-ASCE Joint Committee has broadened — 
this Commentary by including the results of research at other institutions, , 
both here and abroad. The inclusion of this and of contain 


» 


SYNOPSIS 
a _ The evaluation of a considerable amount of res research work has demonstrat- 
ed the applicability of plastic analysis to structural design. . For the type of 2 sa 
BS structure to which its application is intended plastic design results in an over- 


fort open | until December 1, 1959 . Separate discussions should be sub- 


month, a written ‘request must be filed with the Executive Secretary, ASCE, _ Pape wane 
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s __ Among the recent books which should be consulted for a much fuller dis- (i 
cussion than space permits her roblems associated 
ol = stability, connections, and iscussed in Refs.1.1, 
— 


ral steel beams and frames. Theoretical considerations involved in 
the plastic theory and in certain secondary design problems are given. _Ex-— 
perimental verification is provided. . Approximations in the form of 


.. separate and companion report(1. 3) illustrates the procedures : of the 
‘Plastic 1 method of design with ‘specific | reference to building construction, and = 

supplies information to supplement clauses in a specification for plastic 


This progress report covers only a portion of the subject matter. “Other 
papers will follow and will cover “Additional Design Considerations” (shear - 
force, local buckling, lateral buckling, repeated | loading), — at 
41 


Members”, , “Connections” and “Deflections”. 


complete table of contents follows. Progress Report No. 1 includes" 
Chapters 1, 2, 3, a set of definitions, and the Nomenclature for the entire 


a 
— = 


— 


COMMENTARY on, PLASTIC D DESIGN IN ST STEEL 
Chapter 1. 1 Structural Desten 
Chapter 2. BASIC PRINCIPLES ‘Behavior of Material and 


Chapter ANALYSIS. AND DESIGN 1 Assumptions 

3.2 Statical Method of Analysis 


3 Mechanism Method of Analysis. 


‘Chapter _GENERAL PROVISIONS PROVISIONS J Introduction 
4, 2 Types of Contraction 
4. 4 Structural Ductility 
Yield Stress Level 
"Plastic ‘Moment 


Loads 
4. 


5. 2 Beams om 
‘Frames 


Chapter ADDITIONAL DESIGN Shear Force 


CONSIDERATIONS 6.2 Local Buckling 


Wa 20 Reduction of the Plastic 
te to Axial Thrust 
3 Maximum Carrying Capacity of 
Beam-Columns 
5 The Influence of Lateral- 


bs Torsional Buckling 


7.6 Frame — 


Chapter 8. 8.1 Straight Corner Connections 


8.2 Haunched Corner Connections 
on 3 Tapered Haunched Connections 


| 
— 
| 
| 
: 
Chapter 5. VERIFICATIO 
= 


8.5 Beam-To-Column ato Welding 
«8B. 6 Details with Regard to Welding 
Details \ with Regard to Bolting 


— 


9. 2 Deflections in the Elastic Range 
9. in the Plastic Range 
4 Deflection at Ultimate Load 
9. By- Step Calculations 


— 
— 
| 


CHAPTER 1. INTRODUCTION 


eng’ engineering if it can be built it with a 
‘the needed economy and if, throughout its useful life, ‘it carries its intended _ 4 
loads and otherwise performs its intended function. In the process of selecting 


‘The design load of a steel frame may be controlled by a number of criteria, 
any one of which may actually constitute a “Limit of Structural ene i se 


Hypothetical attainment of a specifi yield-point stress (ae- 
2. of maximum plastic strength (assumin idealized behavior) 


| 


- basis for structural design which uses the “allowable stress” concept. Certain 


- provisions also are included in specifications which are intended to insure ‘aaa 


; that the capacity is not otherwise limited. Under certain adverse conditions | ns 


steel may fracture in a brittle manner; and although no specific rule has been 


a available, the - occurrence of such fractures has been reduced to a minimum 


rilabl 
by proper attention to material, design details, and fabrication procedures. -@ 


‘Strictly speaking, a design based on any one of the five criteria tabulated 


above could be referred to as a “Limit Design,” although the term usually has 


_ been applied to determination of maximum strength according to Items 2 —= 
4, _ “PLASTIC DESIGN,” as an aspect of limit design, embraces primarily — 


2 (attainment of maximum plastic strength) and is especially : applicable 
to continuous beams and frames. It is based upon the maximum load aa 
- structure will ‘carry, as determined from an analysis of strength in the plastic | 
4 range (that is, a plastic analysis). Whereas “elastic design” is performed by 


assuming loads and a stress, plastic design is based on 


| 3 


1.2 Plasticity y and Design—Some Advantages and Limitations 


aS <i It has long been known that an indeterminate steel frame has a greater ‘i 
load-carrying capacity than indicated by the allowable stress concept. Such _ 

; frames are able to carry increased loads above the yield value because ~ a le 
pe structural steel has the capacity to yield . While the allowable stress concept 


_ structures has overemphasized the importance of stress rather than a 


is satisfactory for simple structures, its extension to indeterminate steel 4 


asa guide in engineering design. Furthermore it has introduced a com mplexity 


ons is unnecessary for a large number of structures. _ a ae 
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there is no for. an no time the s 


in a steel structure go beyond the elastic range. If this premise were sound, — | fo 
_ then much current practice would have to be abandoned. As a matter of fact 3 ai; a 
it is necessary to consider plasticity in all structural design. An actual fa 
structure is a very complex body with an extremely complicated state of | —  & de 
stress. _ It is an assembly of many individual members joined together to form i m 
a working unit. The individual structural elements such as the beams and .® 
columns come from the mills with residual stresses which are often over one- _ _ 
half the yield stress. connecting the parts together local stresses are -pro- a 
duced by welding and misfits, and there are over- -all assembly stresses. The | ; 4 


structure is pierced by many holes, reinforcements of all kinds are present — 
‘such a as cover plates ar and stiffeners, and many secondary stresses arise owing 
to continuity | , of the structure; for example because of the deformations | 
by the loading, bending and torsion may occur in what are assumed to be 
simple tension members, and axial force and torsion may occur in beams. [2 


& a consequence of these factors (the combination of unknown initial stress, a oa 


ail 


stress concentration and redistribution due to discontinuities of the structure), 
it is inevitable that local plastic flow will take place in the best of designs. . 
An examination of actual load-deflection curves presented later in this repor 

wil demonstrate this conclusively, 
- Numerous examples can be given in which the benefits of plasticity eae a 

: used, consciously or unconsciously, in allowable stress design. By permitting 
a bending stress of 30,000 psi in a round pin, advantage is being taken of the a ‘ 


considerable plastic reserve of strength (70%) that such a shape possesses. _ 7 


ir Similarly the 20% increase in allowable working : stress permitted at points of : 
_ interior support in continuous beams represents another case in which some 

advantage ha has been taken of the redistribution which results from ductility 
_ From the foregoing discussion it is evident that local yielding undoubtedly >, 
occurs in most steel frames the first time full service loads are reached. — 

a Furthermore parts o of these structures will enter appreciably into the ieee = 
before reaching their assumed limit of carrying capacity. 
: ‘This situation will also exist in a structure designed on the basis of plas- we 


ai 4 ticity. Nevertheless it is important to note that at working load, the entesty : 
designed structure is normally in the so-called “elastic” range. A plastic ‘- 
design is based on a load which is equal to the working load multiplied bya = : 
load factor. Any local inelastic deformation at working load occurs only at - 
first loading, precisely as it does in an elastically designed structure. This _ _ 
produces a state of residual stresses such that subsequent | reapplication « of the 
working load | will produce purely elastic stresses. 
e Knowledge of the ductile behavior of steel also permits the designer ae a 
sini costly details intended to provide actual hinges that would be re- as 
_ quired only in case the construction material were brittle. With modern weld- 
- techniques, full capacity splices are often less expensive than mechanical 
hinges or other details that would allow more or less free rotations. Witha — 
knowledge of behavior in the plastic range, the e designer v will realize that the we fe 
elimination of hinges simply results in a stiffer structure with a strength 
- least equal to that calculated for the same structure with hinges. The —— 
in fabrication costs | that result due e to: a better r design of dosreammas based onan 
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ore po exact plastic design which considered strain-hardening would be a 
formidable task. However, if the behavior of structural steel is idealized, * 

go that it is either completely elastic or perfectly plastic, then there is ob- 
tained a simple and yet a reasonably satisfactory approximation termed plastic 
‘design. Within this simplification, the structure is assumed to reach its 

_ “maximum predictable strength at a definite load called the “plastic limit load” 4 
_ or “ultimate load. .” For a large majority of continuous beams and frames, a =. 
design based upon a reasonable factor of safety (load factor) against this ulti- me. 
-mate load provides a more appropriate structure than a design based ased upon ; 
allowable stresses and elastic stress computations. 4 

a _ As was intimated at the outset, there are a number of design factors to ory - 
which consideration must be given over and beyond the selection of members — 

on the basis of plasticity. One such factor is the effect of repeated loading. a 

Another is the onset of buckling before full plasticity is reached and, in 

al, any other effect by which the  Conenee or strengthening effect of defor-_ — 
mation prior to failure is large. fe ‘Brittle e fracture is another problem, of col 4 q 
‘importance both to elastic and to plastic | design. Ultimate load computations 
for continuous frames are based on the assumption that “plastic hinge _ " 
4 moments” are developed at points of f maximum moment i inthe structure and 
maintained during the subsequent loading. design criteria for the stabili- 
ty of details, which under allowable stress design merely guard against the 
initiation of buckling prior to first yielding, require re-examination in plastic © _ 
- design where plastic buckling ‘must be controlled during hinge deformation. a - 

; _ The magnitude of the deflection at working load may constitute still mation, iy 

design criterion. All of these problems are discussed in more later; 


__ A Hungarian, Gabor Kazinczy, first applied the plasticity concepts to the 
design of some apartment-type buildings in 1914.(1.8,1.9) Other early tests’ 
were made in Germany by Maier-Leibnitz.(1.10) Since that time, significant 
contributions have been made to the plastic theory of structures both in this — 
country and abroad. seat. Plastic design is already a part of certain = 

7 cations (it is widely used in Britain) and engineers in this ne: are now 


reghaes all other design procedures. Factors: such as fatigue and buckling 

may become the design criteria. AS an example of a limitation that a oe 

about because of column buckling, most trusses would be excluded from plastic 4 . 
design (with the exception of the Vierendeel type); the method discussed here- a 


4 in requires that “hinges” form, and chord members would not exhibit the 


q are no different from that which exists in any design procedure. ae 


and not the rule. Therefore plastic design is finding considerable application 
in = beams and frames where the members are stressed ions in 


**More precisely | «plastic limi limit design” but referred to hereafter by the 7 


Refs. 1. 4, 1.6 and 1. to he 4%, inclusive 


my 
ress 
t 
ly 
rm 
ro- 
ae 
4 
re), 
ting 
of. 
of 
in 
| 
ic 
yer 
the 
— 
= 
anh 


- _ CHAPTER 2 2. . BASIC PRINCIPLES 


Plastic design takes advantage of an important and unique ecinainy — 


a structural steel, namely its ductility. . Evidence of this’ ductility may be seen | 
by examining a stress-strain curve obtained from a simple tension or com- _ 
pression test. This curve may be represented in idealized form by two ss” 
straight lines as shown in Fig. 2.1. Up to the yield stress level the material 
is elastic. After the yield stress has been reached the strain increases great- 
ly without any further increase of the stress. From this it follows that the _ 
pe attainment of a maximum fiber stress equal to the yield value does not ‘result 
in failure of a bent beam. Rather, the section has a plastic reserve in strength 
2 ped Fig. 2.2 shows the stress distribution at five stages as bending moment is 
_ applied to a member of rectangular cross section. It is assumed that each = 
4 fiber behaves as shown in Fig. 2.1 and that plane sections before bending re- 
‘main plane after bending. The moment-curvature relationship for this beam i 
a is shown in Fig. 2.3; the numbered points correspond to the five stages in cin | 
** Fig. 2.2. Stage 2 corresponds to the yield moment My and stage 5 corre- au 7 
= sponds to the plastic moment, Mp. The exact shape of the moment- ~curvature- 
diagram between stages 2 and 5 depends on the cross-sectional form, but the 
moment rapidly approaches the value of the full plastic moment ae 
to stress distribution 5 (Fig. 2.2). | ‘The * “rapid approach” to this limiting 


is the essence of simple plastic theory. In most 


2. 1. Stress- Strain Curve Uniaxial or 
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‘PLASTIC DESIGN 


Fig 2.2. Successive of Stress Distributio 

the: or curvature relationship is approximated by two straight lines as — 
shown partially dotted in Fig.2.3,.00 
The process of successive yielding of fibers as bending moment isin- e 
creased (stage 2 to stage 5 of Fig. 2.3) is called plastification | of the cross “du 
section. At the section(s) where yielding occurs, relatively large rotations Ro 

a possible without a significant increase or decrease of moments; in other 
words, “plastic hinges” ” develop. The plastic hinge thus formed permits 7 q 
redistribution of moments in statically indeterminate beams and frames. aa : 

Thus, further it increases 5 of the loads are carried by other less heavily stressed 

w parts of the structure, until a sufficient number of plastic hinges are formed 
and the structure starts to behave as a mechanism. Thereafter deflections — 

increase rapidly while the loads" remain practically c constant. 

In summary, a structure will reach its ultimate load as determined 


7 moment has been. reac 


“2.2 Plastic Th ame 


Plastic Theo ry 


In the elastic analysis structure one must consider 


vad 


“Fig. 2. 3. Moment- for Beam in ‘Bending 
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ons _.—s« a'r to form attain the predicted moment and subsequently are able to uacer- Bie 
go sufficiently large rotations. An exception, of course, is the plastic hinge 
which no inelastic rotation is required after the plastic 
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-Continuity—the ¢ deflected shape is assumed to to be a continuous curve, > 
and thus “continuity equations” may be formulated. 

* - Eahieien~th summation of the forces and moments acting on an any 


Limiting Moment—the moment at first yield is the limiting 
In plastic analysis three similar conditions (or : modifications thereof) all 
be considered. regard to © continuity of slope, the situation is just the re- 
theoretically plastic hinges interrupt such continuity, so the 
a quirement is that : ‘sufficient plastic hinges form to allow the structure (or maa 
it) to deform as a mechanism. could be termed a mechanism con- 
i dition. The equilibrium condition is the same as in elastic analysis. Instead 
. of initial yield, the limit of ‘usefulness is the attainment of plastic hinge Ma, = 
moments, not only at one cross section but at each of the sections involved i 
the mechanism motion; this requirement will be termed a plastic moment cone 
dition. A corollary to this is the obvious fact that moments in excess of the 
_ plastic bending strength cannot be resisted. ~The three eras that must 4 
be satisfied | in plastic analysis are, therefore - ial aes 4 


When all three of these conditions are satisfied, it the resulting analysis vee 
for ultimate ‘load is correct because the two limit theorems basic to the plastic 
method are satisfied. satisfied. These theorems will now be d discussed. pal 


general of structures requires consideration o! of the p 
limit theorems of Drucker, | Greenberg, and Prager(1.14,2.1,2. 42.2) whic which have q 
provided c concise and rigorous formalization of intuitively known of 


when changes in are neglected is done in 3 elastic 
solutions) it can shown that failure occurs under constant load and at 
=e stress. In Art. 2.1 this was stated in a different way, namely, that — 

‘the ultimate load is reached when a sufficient number of sections have alee 

their ‘ “limit” or “plastic hinge” moment. nt. Only take 
deformation occurs at constant load. 


theorems are then: | 4 


= or equal to the yield stress, the structure will not fail. At “a 

‘mation, the structure cannot stand up if the rate at which the ex- 


ternal forces do work exceeds the rate of internal dissipation. * 


‘ _ _ Theorem I is a statement of the ability of the material to adjust itself to 
_ carry the applied load if at all possible. It gives lower bounds on, or ae 
values of, the ultimate loading. The maximum lower bound is the plas 
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Theorem 0 is a formal statement of the fact that if a path to failure exists, 


ultimate load. The theorem, therefore, deals with upper bounds on, or — ¥ 
values of, the plastic limit load. The | minimum upper bound is. the 

‘Plane frames and continuous beams are simple enough to permit “ “exact” ai’: 
calculations of the ultimate load. Almost all of this Commentary is devoted sy 


to a consideration of such frames and beams because of the important ‘eaten v4 


they occupy in civil engineering structural practice. 


E 


An actual, therefore three-dimensional, building frame usually is treated i 
s a ‘collection of two-dimensional frames. This procedure is equivalent to  : 
- Setting a number of secondary interaction bending moments and torques equal 7 


zero and is on the safe side according 


_ For. complicated structures it may well be that “exact” | ultimate loads c can- . 4 


_ These two theorems enable the ne of the answer pesytees enough for 


be found due to involved geometry of parts or of the complete ae 


_ The simplicity of plastic analysis opens the way to es plastic hen: 
as contrasted with the trial-and-error procedure which is normally necessary _ 


; in allowable stress design « of indeterminate structures. . Some steps that have 
7 


been taken along these lines will be outlined. 2 —i—i—‘“—s—sS 


_ — There are numerous methods by which the ; maximum strength of a continu-_ 


ous steel structure may be determined. In the . semi-graphical (“statical” or 
-“equilibrium”) method, an equilibrium moment diagram is drawn such that 7 

the moment is nowhere greater than Mp. It thus automatically satisfies the ” 
lower bound theorem. _ The resulting ultimate load is correct only if sufficient 
plastic hinges were assumed to create a mechanism (thus satisfying the upper 7 
bound theorem). ‘In the mechanism method a mechanism is assumed and the 
resulting virtual work equations are solved for ultimate load. This value is Pu 


correct only if the plastic moment condition is also satisfied. = =~ 
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The purpose is to indicate the important of the 


i aaa plastic theory and to describe the essential features of plastic analysis. 
by the “Statical” and the “Mechanism” methods. itis not intended to present ] 
a complete discussion of all the various methods available for analysis and 

design; Refs. 1. 4, 1. 5, 1.6, 3.1, 3. 2 and 3. .3 may be consulted for 


important with ‘regard to the plastic behavior 


of structures according to the “simple plastic theory” are as follows: sine an 


The material is ductile. It has the capacity c of absorbing poate mean 


"mation without the danger of fracture. ah 


~ a moment that is attained through plastic — of the entire ‘cross 
- will occur through a a considerable e angle; in other words, a ‘plastic hinge 


Connections for full continuity will transmit the calculated 


Asa result of the formation of plastic hinges at connections and other 
points where the moments under elastic condition were the critical | 
ones, redistribution o of moment will occur, allowing t the formation of 
oa plastic hinges elsewhere, where the moments under elastic condition 
; 6. The ultimate load may be computed on the basis that a sufficient number 
of plastic hinges have formed to create a mechanism. 
we The deformations are so small that the equilibrium equations can te 


The influence of force shearing force on the plastic moment 


we 


loading is proportional, i.e., the ratios between different loads 


Ln The experimental verification of iain of these assumptions is presented a 
7 later in this report (Chapters 5 and 8). Other assumptions may | need imple- or 


iy mentation to assure that the _ appropriate requirements are met, and this is the 
ioe i concern of portions of Chapter 4 and of Chapters 6 and 7 which will appear in 
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freien in which the number of redundants does not exceed two or three. - - 
& described in Art. 2.2 the objective of this method is to draw a possible 
- moment diagram in such a way that the plastic moment condition is not viata = 
: ed (M< Mp) but with sufficient plastic hinges being formed to create a mecha- < . 
it is desired to determine the total load Wu in span L that may be 
carried by a beam of given moment capacity ina ‘multi- -span structure. One 
of the interior spans A-C is shown in Fig. 3.1la. 
As the load W increases from zero, there are peak moments at the two : 
supports and a peak moment within the span. _In general the moments at the a 
two supports will be unequal, and a typical moment diagram at the elastic _ se 
limit is shown in Fig. 3.1b. The maximum moment within the span is not at . 
= span center line but is near it. As the load increases still further, a . 
greater end moment (C) attains Mp ai and a a plastic hinge forms there. As —_ 
¢ further increases, the other support moment (A) attains Mp, and the Ss 
= moment within the ‘Span i increases as shown by curve 2 in in Fig. 3. lc, 


has shifted to the center of the span and the moment curve 2 is symmetrical. 
‘Finally a plastic hinge also forms at the center, and this corresponds | to the © 


i ultimate load condition. The moment diagram at ultimate load Wu is shown by 


‘maximum ordinate the moment diagram must b be 


a “may be computed 1 from a . knowledge of the yield stress level and the cross — 

geometry. Thus the ultimate load W Wu is ¢ determined. Alternatively, 
‘if the required ultimate load is given by the terms of the problem, the re- = 

ee quired Mp for the beam may be calculated and the beam selected. It will be 

from Fig. 3.1d that the plastic s solution for a ‘uniform, r 


Draw the statical moment curve for load V Wy on a beam of the same . 7 


: ark but with hinged ends. . Draw a new base line for zero moments: ee 
a e”) at the mid-height of this curve. _ Equate the three = 


‘This completes the solution for an interior span of a continuous beam aa a 
restraining moments at least as great as the strength of the beam selected. | 
‘The same reasoning as followed above for a uniform load will lead to the same 
a. solution method for other loading types, such as partial uniform or concen- 


where a smaller beam adjoined an interior span or the case a. 
an end span with exterior not be drawn 
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line reflects the moment capacity at the points of peak moment. In essence, | 
_ then, the failure moment diagram is drawn in such a way that a mechanism _ 
; forms. Curve 3 of Fig. 3.1d satisfies this requirement, of course, because 

the three plastic hinges at a, d, ande result in the mechanism shown in Pe 

a In plastic design of continuous beams and frames it is unnecessary to con- 
sider the elastic distribution of moments such as that shown in Fig. 3.1b or BS i 
the behavior under increasing load. A systematized procedure for the general — 4 
case would simply involve the construction of a determinate moment diagram ra 
(corresponding to curve a-b-c of Fig. 3.1d) which is then combined with a re- 
dundant moment diagram (a-d-e-c) in such a way that a mechanism is formed. ‘. 


at mid- -height of the moment The position of the fixing 


3. 


3.3. Mechanism Method of Analysis 


: Although the matter of personal preference is frequently involved in sli se- 
~ lection of a particular method of analysis, when the number of redundants in a 
frame exceeds two or three the mechanism method usually will be found moet 
suitable. As intimated in Art. 2.2, the objective in the mechanism method is 
to select from all the possible modes of failure the one that corresponds to a 
the lowest possible “ultimate” load. A check is provided through the con-_ 
struction of the moment diagram to make sure that the plastic moment con- 
- dition is not | violated (M< Mp). . In essence, therefore, the mechanism method 
_ gives an upper bound solution that is the correct value of the ultimate load ae 


a Suppose it is required to find the ultimate load he the structure Biss, in 
‘Fig. 3. The section is of uniform moment capacity Mp throughout. 


plastic hinges. These hinges may form at points of peak moment (that i " 
where the shear passes through zero at such places as concentrated load - 


_ points, , connections, ete. ); and therefore the possible plastic hinge locations + 


The next step is to select for investigation the various possible failure 
mechanisms. _ Three of these are shown in Fig. 3.2 . Mechanism 1 corresponds — 
to the action of vertical load P and is called a beam” mechanism. Mecha- 7. 

4 nism 2 corresponds to the action of the norhnceted load P and is often referred = 
_ to asa “sway” or “panel” mechanism. Mechanism 3, on the other hand, is a ~ 
composite mechanism representing the action of both loads; it is a combi- 4 
—s of mechanisms 1 and 2 that eliminates a plastic hinge at section 2. Bea 

3 ect failure mechanism will be the one which results in the lowest a 
since any greater load would mean that the plastic moment condition 
would be violated. The load that corresponds to each mechanism may be com- 

puted by the principle of virtual displacements and is illustrated as follows: * 

_ Referring to Fig. 3.2b, suppose that after the ultimate load is reached the 
= is allowed to move > through : a virtual displacement A. . For pn angle 
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for mechanism 2, 
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0 = 6 My 


Since the lowest is Py, this is s assumed to be the true uitionte load. 
_ It is desirable to check to make sure that some other mechanism was not 
=¢ overlooked. This may be done by examining the plastic moment condition oe 7 
| through construction of the moment diagram. The complete diagram is s shown | 
in Fig. 3.2e, the moment at each section being determined by statics (the 
moments are plotted on the tension side of the member). Tt is found that = 
M< Mp throughout, and thus the answer determined above is verified. It is 
a coincidence that the moment at section 2 is equal to zero. 
_ — The dotted line in Fig. 3.2e traces the moment diagram at the stage at 
ame the hypothetical yield moment would first be reached in the frame— 
| hamely at section 5, where the moment is shown as My. % The first plastic 


hinge would form there, ‘followed by hinges at sections 4, 3, and finally at 1. 


— two previous have in simplified form the most im- 
- portant steps in two methods of plastic analysis. References 1.1 through 1. 6, 
3.1, 3.2 and 3.3 may be consulted for complete apeseennines srucsysepteie 
analysis and design using these methods. __ 
addition to the Statical and Mechanism m methods of analysis, other 
4 niques are available for determining the ultimate load a structure will sup- i 
port. Two of these are the “Method of Inequalities” and the “Plastic Moment | 
- Distribution” method. References 1.13 and 3.4 cover these methods which — 
a particular application for certain structural problems. 
i ~ a Just as in elastic design, where the engineer has available certain charts, 
and with which to analyze standard cases, SO also it 


maf — 


July 
sizes. The designer ‘may use such techniques to shorten the ow even 


o ol but discussion of these aids is not within the scope of this report. 
Many such have been illustrated on design in 
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DESIGN 


— of yield zones have formed to permit the attained whe to deform iia 
ly without further increase in load. It is the largest load a structure will sup- : 
‘port, when perfect plasticity is assumed and when such factors as instability, 


-strain- -hardening and fracture ‘are neglected. 
Allowable Stress Desi n—A design method which defines t the limit of structur-_ 


al usefulness as the load at which a calculated stress equal to the yield point - 
of the material is first attained at any point (usually disregarding local stress 


Plastic Design.- continuous steel beams and which 
"defines the limit of structural usefulness as the “ultimate load.” (The ze 


“plastic” comes from the fact that the ultimate load is Pemgetes Ss from a 
of the strength of steel in the plastic range). wi he 
Factor of Safety.—As used in elastic (allowable stress) design, itisa factor a Gs 
5 by which the yield stress is divided to determine a _ working or allowable stress 
Load Factor. .—In plastic design, a factor by which the working load is ‘multi- 
to determine the ultimate load. choice of terms serves to empha~ y. 
size the reliance ‘upon lo load- capacity of the rather than 


‘ater fiber f first attains yield point — © 


‘Plastic Moment. —The ofa fully - yielded cross-section. 
Plastification.—Gradual penetfation of yield stress from the outer —— " 
towards the centroid of a section under increasing moment. —- i 
“complete when the plastic moment, is attained. 
Plastic Modulus.—The modulus of { resistance to be bending ofa completely yield- 
ed cross section. it is the combined statical moment about the neutral axis J a 
of the cross-sectional areas above and below that axis. 
Shape Factor.—The ratio Mp/ My, Z/S, for a cross- -section, 
- Plastic Hinge. —-A yielded zone which forms in a structural member when the wail 
plastic moment is applied. as if hinged, except that it is 


‘Hinge Angle.—The angle of rotation through which a yielded segment | ofa a beam 


Rotation Capacity.—The angular rotation which a given cross- ~gastionnl shape oF 


can sustain at the panes moment value without prior local failure ~~ 3 
_ Mechanism, — —An articulated system able to. deform w: without a finite i increase in 
load. It is used in the Special sense ‘that the he linkage include real hinges 
- Redistribution of Moment. —A process which results in the successive for- 


3 | EM 3 : 
| — 
| 
| 
of 
of 
— - 
4 
| y 
if 
| 
‘Maton OF plastic minges until tie Ultimate load is Teached. AS a Pesult Ol tie 


arry increased ‘moments. 
pe Loading. —All lo —All loads increase in a ¢ ina ‘constant rat ratio, one to the to the other. 


Area of cross-section. 
€ of shape, Ag= = 
Cross- ~sectional area of plate, = 


= Column height to span n length ra ratio of gable frame (column mw 4 


ie il Distance between « centroids of cross- -sectional area above and ie 


ye Distance from end of cantilever to critical section ie beam. 


» of rectangular cross section 


= Correction factor due to end fixity (restraint) for determining criti-_ 


= Overturning moment parameter (leeward) side. iicastibaaaee 


= Distance between centers of two flanges. 
= Distance two cover 


= = Young modulus of 


Tangent 

Eccentricity. 


= Modulus of elasticity in shear. 


Modulus of elasticity in shear at onset onset of strain-hardening. 


Tangent modulus in shear, 
= Moment ratio in adjacent segment. 


| 
ra 
4 
— 
i, 


= Horizontal reaction. 
4 = Portion of hinge nang that wena in critical (bucking) segment of of 


= Story in multi- story frame. 


= = Moment of inertia; subscripts denote axis. 
= Number of redundants remaining in a structure at ultimate eS 
Moment of inertia of elastic part | of cross-section. 


= Distance from flange face to end of fillet. 

bel 
13508 ‘Plastic moment ratio. 

= = Stiffness factor of a beam. te 


of bar. 


= Leng 
Co = Length of segment (sk (slope- henna: eq 


uation). 


= = Critical length lateral buckling. 
™ ee: = = Critical length (with C = 1.0) of adjacent spans; subscripts | hen and , 


length of connection. 
_ = Length of plastic hinge, 


M = “Bending moment. 
= Critical for lateral of of a beam. 


My at the haunch point. 


moment. 


= Column end moment; useful maximum . 


ssion. 


Moment at load. 


» 
i 
= 
sia 
2 


1988 


of possible plastic hinges. 


Concentrated load. Pe 


= 


Tangent modulus load. 


beam when yield point is reached in flexure. 


 =b/a= roof rise column 


Winer 
Radius of curved haunch. 
Radius of gyration; denote axis. 
Section modulus, I/c. 
“— = Section modulus of elastic part of cross section. - 


Horizontal load applied at eaves which produces overturning 
- ‘moment about the base of mructure equivalent: to that of horizontal — 


distributedload. 
= ee’ thickness; aperene c and t denote compression and 


Stiffener 


External work due to: virtual 
Internal work due to virtual as 


| 


_ = Distributed load per ‘unit of =. 


‘ 


Myc 
3 
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a, 
, 


‘Longitudinal 
Distance to position of plastic hinge ditributed load. 


Ordinate to furthest still-elastic fiber. 
Distance from midheight to neutral axis. oo 
_= Distance from neutral axis al — of half-are 


id 


= Plastic modulus of elastic portion sake cross-section. 


= Plastic modulus of of cross- 


= Trial value of Z, “neglected axial ene. ~ 
= = Lateral 


= Angle t two flanges. 


irtu is men = 


= Deflection. Subscripts, u, w, y denote | work- 


i and yield toad res ectivel rial 


Elongation at fracture 


Measured angl 
ingle change, rotation 

Mechanism angle. 


Poisson’s ratio. 


= ‘Normal stress. 
Lower Yield point. 
= Residual stress. 


= Central angle between points of tangency of curved connection. 

* 

Jé 
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= = Curvature strain- hardening. 


= Curvature corresponding to first ans in n flexure. 
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Report (Proc. Paper 2091) for the Foreword, 
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Discussion open untif December 
Tov. Separate discussions should pe SuD~ 


| 
— ‘This chapter will discuss some of the | basic conditions that Should be si satis- | 
fied in establishing a plastic design procedure. ‘This includes questions re- _ | 
garding types of construction, materials, structural ductility (avoidance of § oa 
brittle fracture), the yield stress level to be used, the plastic moment, the a gl fra 
loads and forces that would be considered as ee to the structure, and the 


In each case the suggested provision will be first, followed by perti- 


st 


2 Types of Construction 


The following types of construction are suitable for plasti 


b) One and two- o-story, , single- and 1 multi-s -span ‘continuous type 

) Multi-story tier buildings with sidesway prevented by walls. — 


- se The necessary continuity may be achieved by welding, riveting, or bolting. ee 


« evices is discussed elsewhere in this report (see Chapter 7). 


on 
a ™~ a Material with the characteristics of ASTM A7 steel for bridges | in 
and buildings should be used, with modifications, when needed, to , 
insure weldability and ductility at lowest expected service tempera- 
itt is s not the intent to specify any one steel, , but to indicate that the important 
property that is required of a material is ductility at service temperature. — is, 
-~ was shown in Ref. 3.1, many of the high strength steels exhibit stress- 7 : mn a 
strain characteristics similar to those of ASTM A7 steel except with a higher : 
yield stress level. It is reasonable to expect that plastic design may be Se he 


BB on to structures in which such steels are used, providing they meet mh 
ahh 
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PLASTIC DESIGN 


i Fabrication processes should be ak as to retain ductility. At | | 


plastic hinge locations unfinished sheared edges and punched holes idee _ 
in tension flanges should not be permitted. Sub-punched and | 
reamed holes for devices would be if ‘ithe 


design, triaxial states of tensile stress up by /geometri- | 


‘This provision toge together with Art. 3. is s intended to assure that 
ian or will not prevent the formation of a plastic hinge. The assumptionof 
the | ie ductility is an equally important aspect of elastic design and numerous design > 
rti- i In plastic design the engineer should be guided by the same principles that | 
govern the proper design of an | all-welded structure designed by the allowable 
: stress methods, since ductility is of equal importance to both. Thus the | ae 
: _ proper material must be specified to meet the appropriate service conditions, 2 
4; _ the fabrication and workmanship must meet high standards, and desi n details 
7 ~ should be such that the material is as free to deform as possible. (4.3) ei oi 


With respect to fabrication, due to the severe cold working involved, Ba, 
_ punched | holes and sheared edges s should not be permitted in in parts that might | 
: be 2 subjected to the yield stress in tension at ultimate load. Punched holes 
_ Would be permitted here if followed by sufficient reaming to remove the cold- 4 
_ worked material. In Ref 4.5 the “effect | of on the brittle 


4.5 Yield Stress Level 
Shear stress, Ty = 19.0 ksi | 


= yield stress level ¢ of 33. 0 ksi corresponds to the minimum yield point 
- ‘permitted in in a mill- -type acceptance test of ASTM A7 steel. Such a test dif- 7 
fers from the test conducted in the laboratory because ofa ‘number of factors, 
one of the most important of v which i is strain rate. An extensive investigation ; 
into the yield stress level has been conducted(4.1 1) using as the test specimen 
| acomplete cross-section of a rolled WF shape. ‘The loading was carried out — 
bs in a manner that simulates “static ” loading. By such a test procedure it was 
‘ possible to include such effects as differences in web and flange strength, 
_ strain rate, and size, since representative cross sections from the very a 
smallest to the largest rolled shapes were included in the program. 
? According to the data available at that time, this investigation 1 showed that 
the most probable value of the yield stress level is 34.1 ksi, with variations 
ranging from 24.6 ksi to 43.0 ksi. (According to the usual are” | 
test, the most probable v value of the yield stress level would be 42.6 ksi). _ 
_ Fig. 4.1 shows the histogram of the ratio of yield stress level 
a stub column test as compared w with the value obtained in a mill- ~type 
_ While 33.0 ksi is the minimum yield stress permitted in acceptance tests, 
mS Aa it turns out that it is very close to the average e basic yield stress: level of this nS 4 
material. the factor of safety includes the 
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minimum. This situation has always existed in and therefore 


4.6 4.6 Plastic Moment Moment 
ay 1) 
= yield stress level 


— 


As pointed out in Art t2. 1, the formation of plastic hinges is of basic im- 
; of _ portance t to plastic design. . Fig. 2.3 shows the characteristic moment-rotation 
- curve of a beam under bending, and the moment at “stage 5” ‘shown in Figs. | 
ey 2.2 and 2.3 is called the plastic moment. It is computed — to ma, <9 
The plastic modulus Z is a geometrical function analogous to ‘sectio 
modulus. It is the modulus of resistance to bending of a completely yielded 
ao _ cross section and is calculated by taking the combined statical moment about _ 
Re e neutral axis of the cross-sectional areas above and below that axis. 
_ As will be evident in Chapter 5, it is frequently observed in tests that the | 
__ moment-deformation behavior is not exactly like that shown in Fig. 2.3 (See 4 


* Because of strain-hardening, the resisting moment is 


= 


Fig. 4. 1. Frequency D Distribution of the Ratio (Stub- -Column) 


EM 
sin 
| 
ig. 9.4, lor example). 
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greater than the value to Eq. ad. no 
_ simple theory can take this additional reserve in strength into account without | 


_ For material whose characteristics are not similar to A7 steel, but gall 
exhibits continuous strain-hardening, it might be desirable to arrive ata iar. 


_ semi-empirical value for the “plastic hinge” moment. Studies would have to ‘* 


be made on the particular material (including bending tests and tests of inde- 
terminate structures) to arrive at for the plastic 

& _ The loads to be provided for (allowable loads) should be those — ee 


that are customary for the particular type of construction. | These 


loads are increased by a load factor to obtain an ultimate load- ua 
ing. Members are selected on the basis of their plastic strength 


| 

A margin of safety is in elastic design by the use of allowable 


tba unit stresses obtained from a unit stress level assumed to represent failure = 
4 da g which has been reduced by a “factor of safety.” In plastic design safety — cS 


j > is achieved by multiplying the given service loads bya a “load factor” as dis- a 
| The use of plastic design does not involve any changes in the magnitude 
a the service | loading De specified for a given structure. The difference is that, 
“1, 3 the case of plastic design, members are selected so that the structure will ne 
support the computed ultimate load P,. In elastic design the members 
— are so selected that allowable unit stresses ' will | not be exceeded at ‘service — iv 
phe =~ Py . (Asu used here P is the critical combination of given aaa” = - 4 


' variable loads used as the basis for the design by either method.) .— ae 
‘The conditions that would be investigated for building construction 


2. Dead load, plus live ‘load, plus wind or earthquake on, en 


il a is assumed that the live loading is static and proportional, ies for the 
_ characteristic fluctuations of live load found in buildings. ae unusual con- — 


ditions, deflection | stability would be investigated (see Art. 6.4) 


a. In the design of structures to resist the dynamic loading resulting from —— 


J 


blast forces, plastic desirn concepts generally are u used. be However selection 
the loading and the typ. of design are largely a matter of judgment and are 
| _ based on studies of vulnerability, consequences of failure, and required radi- i . 

ation | protection. In making the design it is necessary to consider the nature =e: i 


structure, as a as factors not considered in 


— 
ation 
| 
_ 
if 
| 
“ah 
| 
— 
UE UIE mcrease am eld strengu Ol tne material, ElleClive 


| 


‘ equivalent lateral static loads to produce a desired level of earthquake re- 
‘sistance. In conformity with such practice, the specified loads can | be | multi-— . 
plied by a load factor to yield design ultimate loads. 
5 _ Concepts of plastic analysis are currently being considered for design to 
‘resist earthquake forces using procedures that take into account the —_ 
and random nature of earthquake loading require that consideration be given 
to loading approximations used in the design, as well as to many of the factors 
noted above as of importance in designing for blast forces. Greater ne 
_ finement in design procedures using these concepts must await further eet. 
‘The complex nature of such an approach makes it almost mandatory to us Ne 


digital com uters in the anal sis. 
[Dead load plus live load, wi 


— 


"sign and the load factor used in plastic design are not concerned d alone with 
the possibility of overloading. _ Other : factors which influence the selection of 


Approximations and uncertainties in the method of analysis 
Presence of residual stresses and stress concentrations* 
Under- -run in physical properties of material nee 
& 5. Under- run of cross-sectional dimensions of of members 


Lo Location and intended use of 
upon the) type of and its intended use, the importance 


7 of any one of the above factors, as compared with the others, can vary some-_ 
; 3 what. One might arrive at a precise over-all load factor F in each case hae 
; sufficient statistical data were available to weigh properly the importance of V 
its various constituent parts; but any resulting departure from current ae 
: practice would be equally applicable to the allowable stresses specified for 
_ use in elastic designs. Since such statistical data are not available, it would | 
seem logical to draw upon the vast experience gained from allowable _— 
to obtain a single value, the range of 


‘ portioned, there being no necessity for requiring any greater margin of cntety. 
_ merely because the structure is redundant. | For simple beams the load factor | 
iE equal to the ratio of the ultimate load Py to the working load Py; thus F 
equals P,,/Py. Since here the moment varies linearly with the load 


at 


*While this factor contetbutes to as to the | recise stress level, 

a the discussion in Arts. 1.2, 5.1 and 5.3 shows that it may not influence _ 


a 
ite 
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re fis the shape factor. 


_ Since adoption of its Standard Specification for Structural Steel for Buildings 


tent even before that, the basic allowable working stress in building design in 
_ the United States has been 0. 6 of the specified minimum unit yield stress of | ca 
the steel furnished. Restated, the load factor —_—" the —— 
ultimate capacity of these beams has not exceeded 


30,000 psi 33,000 psi 1.65f 

ci “ome The formulation of a satisfactory load factor is therefore mie only 

“upon the determination of a shape factor representative of the simple andl 
‘The variation of the shape factor for WF beams and columns, andfor 
American Standard beams is shown in Fig. 4.2. For WF shapes normally | used 
_ as beams listed in the “section economy” ‘table(4 -4) the shape factor varies 
from 1.10 to 1.18 with an average value of 1.134 and a mode (most frequently 
_— cbeerved value) of 1.12. For WF shapes normally used as columns that ap- 

_ pear in the “column” tables of Ref. 4.4 the shape factor varies from 1.10 to Ry 
x .23 with an average value of 1. 137 and a mode of 1.115. The shape factor 


Fig. 4.2. 
8. 


pending on the choice of the shape factor. 


SHAPE FACTOR FACTOR OF SAFETY* 


alue ane 
— Mode WF beams 


Columns 


* 18 — Average for American ce 


I beams = 


‘The two most reasonable values for the load factor are 1.85 and 1.88. The 
former would seem appropriate because it represents the shape factor = 
recur most frequently i in beams. value 1. 88 implies an in 
_ In the case of gravity loading in combination with wind or earthquake - 
forces, stress design specifications permit a one-third 
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. The Variation of Shape ‘Factor for WF Beams and amare F 
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computed stresses. - Consistent v with | this ; allowance, the value of F for co com- 
dead, live, and wind or would be 3/4 x1. = = = 1. 40. 
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CHAPTER 5. _ VERIFICATION OF PLASTIC THEORY 


+ 4 In ‘this chapter it is shown that hina actual behavior of structures under test 
verifies the predictions of plastic theory. In Art. 5.1 it is demonstrated that | 
structural steel exhibits the ductility assumed, and that plastic hinges will ee 


hem and allow the ‘necessary redistribution of moment. Article 5.2 presents 
the results of continuous beam tests, ts, and shstainiasnindiaii Art. 5. 3 shows ho shows how tests of 


q 


Fig. 5.1 shows the tensile stress-strain curves obtained from. two wo coupons 


cut from two separate locations of an 8WF40 beam. They are typical of the _ 
ensued of ASTM-A7 | steel. The steel deforms recears about 15 times the 


a 
— 
— 


ings 
strain at the elastic limit and harden. * Although the 
s) 5 data is plotted well into the strain-hardening range, the strains shown are ~~ 
t still considerably less than those at ultimate strength (tensile strength). The 7 
compressive and the tensile stress-strain relationships are quite similar. -_ 
bh fact the properties in compression are practically identical with those in 
tension. Fig. 5.2 shows in idealized form the stress-strain relationships tor 
ASTM- -AT steel drawn according to average values obtained in 
The Plastic Moment and the Plastic Hinge 
oe As a demonstration that the plastic | moment is attained through einen 
Cation of the cross section, Fig. 5.3a shows a 
a beam, a ‘Portion of which is in pure bending. ( The dotted line is the q 
‘idealized curve and the solid line through the circles shows the results ofa " 
test. The theoretical stress distributions (according to the ‘simple 
a theory) at different stages of bending are shown in Fig. 3b. Below these in 3 
‘Fig. 5.3c are shown the corresponding stress distributions as determined 
from SR-4 gage measurements. It will be seen that plastification of the cross 
_ section does occur, and that the bending moment corresponding to ‘this con- - 
dition is the full plastic ‘moment as computed from the equation Mp Bey. 
Although there will be inevitable minor variations from the result 
in Fig. 5.3a the many tests conducted on rolled indicate that most | 


*ASTM AT requires an elongation in 2 in. n. of ‘not less than 24% at failure, ‘ade 


(qRical M-¢ curve obtained from 


= 
is 
7 Be 
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a To be sure, a somewhat unrealistic loading condition has been taken since’ = 
“pure moment” is a condition not likely to be encountered in actual structures. - 
7 Usually there will be a gradient in moment, as when a single concentrated load _ 
me is applied to a beam. In such a case the deformation tends to be concentrated - 


= the load point (the point of maximum moment). e Because the plastic on 


deformation is more localized, the strain-hardening region is reached at a 
lesser deflection; consequently, the beam tends to develop a moment greater bs 
than the plastic moment. Typical of the behavior of a beam under moment | 
_ gradient is that shown in Fig. 5.4. (5.2) The theoretical load-deflection curve 
_ with and without the inclusion of strain-hardening is shown by dashed lines. 
The results of a test of such a beam are shown as a solid line. it will be noted — 
that asa ‘result of the strain- -hardening phenomenon, there increase in 
load-carrying capacity as the deformation continues beyond the yield level 7" ; 
The decrease in measured resistance, which occurred after a large hinge ro-— a 


tation at the center, resulted from local eas of the flanges - followed by 


Fig. ; 5.3. Experimental Verification of the Moment-Curvature Relationship 
_ Theoretical anc Computed Stress 
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5 ‘PLASTIC DESIGN 

‘Thus. -hardening improves the moment- carrying capacity of a beam. 

- - Although it is neglected in the simple plastic theory (except for checking a 7 
beam for stability against buckling) this additional reserve strength is 

| present in most ordinary structures, and this contributes to an actual — 


of safety that is the value assumed in the sim simple plastic 


a; upon to form at connections and at concentrated load points. a Development | of 


the plastic moment is one of the sources of reserve strength in structural 
wut steel beyond the elastic limit. Another source is the redistribution of moment 7 
tin Fig. 5.5 is shown a picture of the e redistribution pr process— —as predicted 


beam to simulate the condition of third-point loading on a fixed- ended —- ; 


and as obtained experimentally. A test was made on a continuous 


ati span;(5. 3) thus experimental data was available to compare re with the theoreti-_ 

cal predictions. fixed- ended beam and its various are 

age 1—at the computed elastic limit 


after the plastic hinges have formed at the ends’ and the load has 
_ increased towards its ultimate value 
heal 3—when the theoretical ultimate load is first reached, and — or a 
Stage 4—after deformation has been continued through an arbitrary ad 


The figure shows (a) the loading, (b) the deflected : shape at the four —_ 7 


near the wap (e) and at the | center (f). 
_ In the elastic range (stage 1) it will be seen that the beam behaves just as as 
we ‘assumed by the theory, the moment at the center being one- -half the moment 
_ at the fixed ends. . (Figs. 5. 5c, 5. 5e, 5.5f). AS the moment at the ends ap- 

4 


7 : proaches the yield moment, the curvature ¢ commences to increase more 
me - Hapidy a plastic hinge begins to form (Fig. 5.5e). Because of this “hinge 
oe action”, the additional moments due to increase in load are distributed be- 
— the ends and the center in a different ratio beyond the elastic range than 
7 before. As long as the beam is is elastic the increase in moment at the center 
_ corresponding to a load increment is one-half the increase at the ends. How- 
ever, after a plastic hinge forms at the ends, most of the increase of moment _ ts 
— at the center; the moment increment at the ends is “el (Fig. 5. 5e, Dy 


5. 5f) This is the process known as redistribution of moment. a 


what more flexibly than before. (Fig. 5.5d) “At stage 2 the elastic moment | 


capacity near the center is practically exhausted. ‘It is quite evident from 
_ Fig. 5.5 that substantially all of the moment capacity | has been absorbed by the s 
time stage 3 is reached (ultimate load). Beyond this, the beam simply de- a 


plastic hinges at the ends and center rotating further, 


am Clear evidence is therefore available that redistribution | of moment occurs 


— 

| 
ug ation Dia nges, allowing the structure to reac and 
ultimate load predicted by simple plastic 
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Incidentally, y, sketch of Fig. 5 the from 
_ the elastic to the > inelastic range that is typical for continuous steel beams _ 
a and frames. . Theoretically, upon first loading the structure should remain _ 
elastic up to stage 1. However, because of the combination of unknown initial | 
‘stress conditions and discontinuities, local plastic flow takes place at a lower 
bh load than that which . corresponds to stage 1. But there is no effect gece 


fabricated from rolled sections. structure and loading are shown 
. scale at the left. Next, the size of member (or members) is indicated. To- ar 


graph is ; the ‘computed elastic limit (on a non-dimensional basis) and the « end 
of the shaded portion is the observed maximum strength, 


Particularly remarkable the continuous beam tests of Fig. 5.6 is one 
10 


conducted by Maier- Leibnitz (1.1 shown as the next to the last structure. he 
this experiment, prior to applying the vertical load, he raised the center sup- 
port until the allowable working stress was just reached, with the result 1 that — 
: application of the first increment of external load was, in fact, a load greater 
_ than allowed by the specifications. In spite of this, the computed ultimate load 
was attained. The observed ultimate load in this test was 3% of that 
a of the two structures shown immediately above in Fig. 5.6. Adega : 
A The continuous beams shown in Fig. 5.7 were tested to show that members 
of otherwise inadequate strength may be en achieve th the — 


this country as part of the experimental verification of the plastic theory. “gt 
The span is 40 ft and the frame was fabricated of 12WF36 rolled shapes. Not 
only has the computed ultimate load been reached at the stage shown inthe _ 
photograph, but the frame has absorbed considerable additional plastic defor-— 


= mation while sustaining a load slightly , in excess of Py. Fig. 5.9 shows the 


_ load-deflection curve for this structure. The dashed line shows the predicted 


behavior based on theory and the series of open circles connected by the solid 


_ line represents the observed behavior. This result demonstrates once again 
- that inelastic action (due to local effects) commences at a load considerably - 


ol less than the predicted yield value. It also shows that the ultimate | load is not 


va It is of interest that at ultimate load the excess 0 of ‘actual deflection shove jj 


Pe ‘a ‘a computed value for this frame was no greater than that observed at the ¥ 


predicted yield load. means that the methods for computing such | de- 


flections are as dependable as the elastic deflection calculations. 
a Figs. 5.10, 5.11, 5.12, and 5.13 show frames tested both in this country and 


a | um a ap represent some of the structures which have been tested to ices 
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en end of that part being the observed load. a agreement is 
§ served at maximum load except for those cases in which strain-hardening ac- 

counted for an increase; it is better, in fact, than the agreement between the — 
predicted load at first yield and the observed value. In Fig. 5.10 testing of the 
fourth frame was interrupted in order that the fifth test might be carried out 
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the same structure but with a different of ‘horizontal to vertical 
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‘Fig. 5.12. 
wi Predictions of Plastic Theory 
I 


n view of the notable ‘agreement between plastic ‘theory ar and the results of 


these tests, the applicability of the plastic method to structural design = 


involving continuous steel beams and frames is is 


Respectfully submit ‘submitted, 
A. Frank Baron 
J. M. Crowley Daniel C. Drucker 
Samuel Epstein 4.2 
LaMotte ¢ R. Higgins — 
William H. Jameson Harry N. 
Bruce G. Johnston Eivind Hognestad 
Bruce G. 
William H. Munse, Jr. 
Egor P. Popov 
= George Winter 
R. Higgins, Chairman Lynn S. Beedle, Chairman 


William J. Hall 


Steel Engineering Mechanics 


Welding Research ‘Council Society of Civil Engineers 


Lehigh Project Subcommittee Committee on Plasticity Related to 
Le 


STRUCTURE, LOADING, AND SHA 

a 

aie 
»f the 
| | 
— 


rr. 
g0___100 


¥ (two frames y 
tested 
poraliel) 


= 


q 


frames 
tested in 
poraitel) wall 1 


Big. 8. 13 of Frame Test Results with 


— 
| 


PLASTIC DESIGN 


Beedle, Huber, A. W., Stress and the Compressive 
> _ Properties of Steel”, A Summary Report, Fritz Laboratory Report Lae 


4,2, Huber, W., “Residual Stress Measurement”, Fritz Laboratory 


220A. 17, (March, 1955). — 

Thurlimann, B., “Modifications to the Plastic Theory” , Proceedings of f 


‘the AISC National Engineeri Conference p. 50, (1956). 
AISC, “Steel Construction”, American Institute of St Steel. Construction, 
New York, 5th Edition, (1950). — 


4 5. 7 Harris, L. A., Newmark, N. M., “The Effect of iain Edge Con- | 


ditions on Brittle Fracture of Structural Steel”, A.R.E.A. Bulletin 59 
(538), | pp. 247- 289, , (September- October, 1957). 


ding Journal, 27(11) p 


Driscoll, G. B Beedle, I L. 


Yang, c. H. S., , Johnston, B. “Residual Stress an and 

Compressive Strength of Steel Beams”, Welding Journal, 31 1 (4), p. 


Hall, , “Shear ‘Deflection of Wide- -Flange_ 
Beams in the Plastic Range”, Transactions of the ASCE, , 122 22 p. 666 Ps. ie 


the 

fs: is 5. 5. Popov, gE. Pp, Willis, J. A., “Plastic Design of Cover-Plated Continuous 4 
a Beams”, ASCE Journal, 84(EM1), Paper No. 1494 (January, 1958). ee 

., Johnston, E. R., Beedle, L. 8 Ss. Welded 
Portal Frames Tested to 


5. “a 
Single-Span Frames Under Combined Loading”, 


Driscoll, G. C., “Test of Two-Sp Gabled Portal Frame’, AISC 
National E ineeri Conference Proceedings. 74 (1956). 


: - 9. . Nelson, H. M., Wright, D. T., Dolphin, J. W., “Demonstrations of Plastic 
Behavior of Steel Frames”, “ASCE Proceedings, No. 


10. . Baker, J. F., Roderick, J. W., “Tests s on Full- ‘Scale 


_ Proceedi s of the Inst. of Civil Engineers (January, , 1952). ae 


| 
| 
= 
— | 
— 


July, 1959 
.11. Baker, J. F., Eickhoff, K. G., “The ofS Tooth 
Frames” Conference on the Correlation between Calculated and 
served Stresses and in Inst. of Civil Engi- = 
p. 107, (1955). ai 

5 12. K. = “A ona Pitched Roof Portal”, 
Publication, IABSE, Fi ifth Congress, Lisbon 


= 
7 


eT? 


| 

i= 


oOo 

_& 


= 
— 


July, 


“an 


at 


ENGINEERING MECHANICS DIVISION 
i Proceedings of the American Society of Civil E ngineers 


STRATIFIED ‘D FLOW IN INTO SINK 


a fluid (with a stable linear sity variation nat 


ally investigated. The experiments show ala when the densimetric Froude 

number is less than 0.28 the flow pattern is divided into two horizontal : 

regions: an upper, essentially stagnant, region and a lower region in te 

_ the entire discharge is concentrated. _ The specific results obtained are useful 

for the solution of many engineering problems stratified flows. = 


INTRODUCTION 


_ The experimental investigation ofa stratified fluid flowing along a hori- 


zontal channel into a line sink was undertaken to gain additional insight into 
_ the mechanics of stratified fluid flows. The practical application of an 
= problem, and those similar to it, is encountered in the cooling of thermo _ ; 
fl electric generating stations. . In such installations it has been observed: i 
the coolant discharge water flowed, by virtue of its lower density, into the | 
proximity of the coolant inlet where it was subsequently pumped through & 
cooling system with an attendant loss in operating efficiency. 
2 _ Other examples in which stratification phenomena are present are in the o 
drawing off of crude petroleum from underground reservoirs, and the removal 
of salt water that has encroached upon a supply of fresh water. ‘This latter 
ae has assumed considerable importance to agriculture in Holland. _ me 


analysis of “Yin(1) for the ‘type of flow under consideration the 


that there is no tendency for flow (i. e. , discharge of ‘the 

Note: Discussion open until December LZ 1959, To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 2093 is 4 
a part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
- the American Society of Civil Engineers, Vol. 85, No. EM 3, July, 1959. aa) a 


Prof., Dept. o of Eng. of Ann Arbor, 
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for Froude numbers equal to or less than a critical value of 1/7. His results 
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"heavier fluid to the exclusion of the lighter for Froude 
greater thai than 1/m. Any separation would be expected to occur when the seus 
2 number was equal to or less than this critical value. ‘Sore: aie 
‘The present work was undertaken to determine the critical Frente number 
and to compare its mi magnitude wit with the value obtained 


— The mathematical solution for the two-dimensional problem em of a stratified 
a fluid in a horizontal channel flowing into a line sink was obtained by Yih. (1) 
An abridgement of this work follows in order to give a complete background eit 
for the experimental work that was undertaken. 


“. al The density variation in the channel at a section far from the ‘sink is con- 


Lae ‘sidered to be linear and of the form: 


_ ‘in iil y is measured vertically at the bottom, do is the depth of fluid in 7 
_ the channel, po i is the density at the bottom, and pg is the density at the top. _ 
. Ih this discussion x will be measured horizontally from the sink and only Paar 
_ steady flows are considered. The fluid is incompressible and changes in > a 7 
due to diffusion are excluded as being negligible. 


49% 


in which u and v a v are the respective velocity components in the x and y di- Pe 


rection, and the subscripts denote partial differentiation. the 


 contionity equation to be written inthe usual form 


and the use of the Senction W so that 


"isintroduced suchthat = 4 


7 ‘The e uations of motion can then | ~ saree in terms o ‘ 


tangular cartesian coordinates, 


is the acceleration constant, and Hy ( isa that must 


| 
in 
vas separation with the Froude number, and also to observe the flow patterns 
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STRATIFIED | FLOW 


This is accomplished by utilizing the conditions that exist 
far upstream. If the fluid originates from a large quiescent reservoir and rs 
flows s horizontally into t the channel, , it can be owe that far upstream 


in which the subscript ‘designates the the andAisa 
constant. One can also deduce if is the far 
upstream in the actual flow, then 


« 
— 


and 


‘This follows from the fact that the solution yields an edie); shove the line 
sink, which becomes more and more elongated as the Froude number is re oa 
duced. At F = 1/m the eddy reaches infinitely far upstream and invalidates _ “ 
the boundary conditions necessary for the evaluation of H;(¥1) in Eq. (6). The 
flow pattern for F=.35isshownin Fig.l. cr, 
‘The conclusion is reached in the ® original analysis that when the Froude ih 


i. flows at Froude numbers > anne or less than cr remained to be answered. 
this question an in experimental program initiated. 
Analysis 


‘the case of stratified flow for which there is a of 


pattern into two regions—an upper, , essentially stagnan 
region in which the flow is concentrated—the following variables can ~ oe 


selected oblem (see Fig. 9): 


oll 
hat 
— 
ed | 
This equation along with the appropriate boundary conditions can be solved 
(10) 
(4) 
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=. 


5 CO 


= ‘Depth of liquid in the ne channel ad 


= = Density of at the channel 
a = Po - P2= Density difference between the liquid at on ines: of the 


"= 45 
- Gravitational constant 


= Viscosity « of liquid in n channel ute ait 
one term included s since for small Ap the differences 
in viscosity throughout the liquid are exceedingly small and can be neglected. , 


in which $1 i is some me function t to t be determined. By in itll: the der densimetric 


— |“ 
Streamline pattern predicted by analysis for F * 
q Volumetric discharge per unit ft 
= = Heicht of zone in which flow is concen eq 
| 
«CS 
a 
a! 
: 


The ‘experiment shows Fo to variable in the function. = 


_ parameter Ap/po is important only when it is not small, and the Reynolds — 
number, GP,/ is of secondary importance for the type of flow under consider- 


‘¢ In the section devoted to a theoretical development and where separated 
flows were not considered, oe: Froude number —_ defined by Eq. (12) as 


- governing the flow. If the derived velocity distribution given by Eq. (8) is used 


and small differences in density are considered, the | constant A can be elimi- 


it to as the characteristic length the height of the 
_ region in which the flow is concentrated, , dy, rather than the channel | depth, ma 


results, in B is upstream density gradient defined by Eq. (1). Fy 
' related to Fo by the multiplier (dy/dg)2, and the significance of considering | 
Fy as a variable will be seen when the aa icceait results are discussed aa 


5 _ The experimental work was carried out in a flume at the Hydraulics Labo- at 
ratory of the National Bureau of Standards. This flume is 18 feet long, 2 feet 
high, and 1 foot wide (see Fig. 6). An eight foot long partition was placed in — 
the flume in order to obtain a uniform channel of 6 inch width for a consider- _ 

able distance upstream of the horizontal orifice. In this } Way it was possible 
to reduce the total quantity of discharge and to minimize the effects of the — 

associated decrease in head pressure. The discharge orifice (see Fig. 6) with 

its rounded entrance | to eliminate the formation of eddys was 11/16" in width j 


i 


—— 


and is s judged to be narrow enough to approximate the line sink considered in 

a The flume was filled with 15 distinct layers of salt-water, each with a —s 

different density, to a height of 22.85 inches (58 cm.). The alternate layers of 

water were dyed with a small amount of nigrosine prior to admission into the > 

_ flume so as to delineate the flow pattern during the discharge of the flume. Ne me 

| This technique for marking the flow was very satisfactory, and was prompted 

_ by the theoretical consideration that the streamlines would be lines of al 

density. The densest solution was first added to the channel and then each ~ 

lighter was onto the previous one, 


| 
-€ |] which is the same result as that obtained by dimensional analysis consider- a 
is 
4 
— 
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Fig. 2. Flow pattern obtained for F= 
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Fig. Photograph of flume immediately pr prior or to in aye in the 
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5. Flow p pattern for Fj = 0. 245 (Test 5). 
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used for this en process is ; shown in Fig. 1. win this mente: there ere 
resulted a stream of liquid that was uniform across the breadth of the flume 
and which, , when the rate of liquid addition v was sufficiently low, yielded alm almost 
A period of about 18 hours elapsed between the time that the flume was 
completely filled and the time that the experiment (i.e. emptying the flume) 
was performed. This period of time permitted the diffusion process to smooth — 
out the abrupt density differences between adjacent layers and | Shaseky promot- 
ed the > establishment of the desired linear density distribution. . To this end 
the layers were made as narrow as practicable, 4 cm. . However, the top and 
bottom layers were made 3 cm. high in order to compensate for the absence > 


of diffusion at the upper and lower surfaces ofthe channel. =” 


a After the overall density difference fora particular test was decided upon, 
the prescribed density of the individual layers was obtained by adding quanti- 


7 3 ties of sodium chloride, the required weight of which was determined from 
: appropriate tables. An effort was made to obtain a qualitative picture of the 7 
= density distribution after diffusion had taken place by using an electric probe | a 
that measured the conductivity of the solution between two platinum electrodes. 
2 3 & This apparatus | was calibrated with solutions of known salinity and the current — 
recorded for a given voltage. A traverse was then made of the filled flume to 
determine the density profile. The results obtained with this probe verify 
the establishment of a continuous linear density profile from the original 
"discrete or staircase distribution. Fig. 8 the results ofa 
sity traverse made before one of the experiments. The change in the 
distribution —from. staircase to linear—is clearly | evident even though the 
/ = second set of density (i.e., conductivity) readings is consistently greater than ] 
gq &§ the first set. The cause for the apparent, but unrealistic, increase in density — an 
sj 2 between the two traverses was investigated but no complete explanation was Bi 
% g obtained. With the linearity of the density profile established, the magnitude as 
8 of the gradient used in the computations was determined from the channel - a ; 
= depth and the difference in density of the liquid at the top and bottom of the | "i 
ae g channel. The values for these two densities were obtained by knowing the it & 
 § weight of dissolved sodium chloride in the top and bottom layers a and using the zi 
aforementioned tables which relate solution density and the amount of solute 


_ © § The flume was emptied by opening a gate valve that was located down- Cie 


ba 2) stream of the sink. * Different valve settings | gave the various s discharge rates 
a used in the experiment. The discharge rate itself was determined by observ- — 
ing the time required for the free surface in the flume to drop ic 1 cm. in in height. 
its ion 
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STRATIFIED F FLOW 


= iF The liquid mons was measured with a multi- -pronged hook g gauge, e, and a record- 
Ang ing chronograph was used to note the various time intervals ee 4 
= _ Photographs of the flow pattern were taken in the vicinity of the discharge — 
all orifice whenever the free surface dropped 1 cm. in elevation (see Figs. 2, — 
A ‘ 4 and 5), and motion pictures were taken near the upstream end of the flume - q 
a in order to record the flow “at infinity”, 
Ce The alternately dyed layers were effective in showing | the flow pattern but 
: they did not give much information regarding the magnitude and direction of — 
_ the velocity components within the layers. In order to overcome this diffi- 
“culty, vertical dye streaks were introduced by dropping potassium permanga- 
d nate crystals into the flume just prior to the experiment. These dye trails © 
along with a grid placed upon the side of the flume enabled the experimenter 
to ‘observe more thoroughly the motion | of the liquid. 


be The most significant venite ene the experiment was the discovery of the fi; 


separation of the flow pattern in into two regions for densimetric Froude 
"numbers below : a critical value. division in the flow pattern was hori- 
~ zontal and extended the full length of the flume. The upper region contained — 

liquid which did not flow out of the discharge orifice > while the lower region :. 

‘consisted entirely of liquid that flowed out of the channel. As the ony 
"number was decreased the region occupied by the upper, or relatively static 
portion became increasingly greater. . The pertinent data for the six tests that 
were conducted are contained in Table I and summarized in Fig. 9. Sade) 
‘The location of the dividing streamline was obtained by examining and 
¥ measuring the photographs that had been taken. Several successive photo- 
. graphs for each test were used which accounts for the presence of more than — 
a one data point in Fig. 9 for a particular test. The data presented were gotten 
es from | those photographs that were taken after the unavoidable inertial effects — 


SUMMARY OF EXPERIMENTAL DATA 
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in the neighborhood of ar However, the quantity of flow that was needed to — 
achieve this Froude number was so great, and the accompanying drop in head 
so fast, that an extensive investigation was not considered justified with the 
- available —— > Therefore, only a preliminary experiment was made for 


Yih for F = "35 (Fig. 1) and the experimental result for F = .347 (Fig. 2) shows | 
a striking agreement. The bulk velocity is 0.31 Lp. s. in this photograph which 
was taken 35 seconds after the initiation of flow, 
_ The data presented in Fig. 9 indicates that separation of the flow pattern | 


_ commences for a Froude number of 0. 28. This agrees quite well withthe _ 


Ys These components were small in magnitude and would reverse their direction 
as the flume was emptying. ‘The direction of the motion was occasionally a 4 


% _ rameter when dealing with stratified flows of the type under discussion. 


4 _ large to conclude anything significant regarding the velocity law. The motion 


value of 0.318 (i.e., 1/m) suggested by the analysis as being the critical Froude 
_humber—in the sense that a significant change of flow pattern is to be meee 
ed. The experimental curve has been drawn to pass through the origin. “This: } 
= be expected, not only from the data, but also from another consideration. 
Yih(2) has shown analytically and experimentally that as a horizontal dis- 
turbance in a stratified fluid becomes infinitesmal in strength, the vertical ex- 
- tent of the induced fluid motion is confined to a region equal to the vertical = 
height of the disturbance. | Thus, as the discharge is reduced to zero, it can be 
c expected that the region of flow, hence > dy, will approach the e height : of the line 
Fig. 9 also shows the curve (dy /dg)2 = mF. It is proposed that this curve 
_ represents the results that would be obtained with an inviscid fluid. The two 
_ end points of the curve can be justified, and the form of the curve is consistent 
a with the observed data and the expected influence of viscosity on the flow. The 
_ proposed relationship also supports the conjecture that Fy would be nearly 


constant if plotted against dy dy and that Fy may be the more ‘significant pa- 


_ Fig. 10 represents the accumulated data with Fy as the independent variable. by 
a It was observed in the course of the experiment that there were horizontal 
velocity components in the upper region (i.e., above the dividing streamline). 


opposite at different elevations in this region of the flow. There was some 
indication of wave motion within the upper portion, but there was no evidence | 
_ whatsoever of mixing of the liquid. _ Thus, compared to the lower region, the 
upper is considered to be 
_ With the motion pictures it was intended to substantiate the upstream ve- 
locity distribution utilized by Long(3,4) and Yih in their theoretical analyses. 
However, the differences in densities that were employed were not sufficiently 


pictures did verify the region of flow separation as shown by the still sata 
taken the other endofthe flume. 


CONCLUSION 


ie It has been aideaieaas for the flow of a stratified fluid with a linear nes 


D : variation into a line sink at the end of a horizontal ‘channel, that a 


critical densimetric Froude number exists, below which the flow pattern di- 
vides into t two distinct regions. | There i is an upper region composed of fluid = 
does not flow into the sink which remains relatively 
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a all of the fuid is As 
_ Froude number is decreased to zero the height of this region of discharge, dy, 
is also reduced to zero. The variation of d, with the Froude number is given _ 

; in Figs. 9 and 10 which show the critical Froude number to be approximately 


a ey The experimentally determined critical Froude number is somewhat less _ 
= € than the theoretical value of 1/7 based upon the assumption of an inviscid _ 

fluid, but it is sufficiently close to give support to the eoretical analysis de- 
velo) ed for higher Froude numbers. = 

4 The experiment on which this paper is based was by the ong 
the Hydraulics Laboratory, National Bureau of Standards, Washington, Cc. 


Without the kind permission of Drs. G. Schubauer and G. Keulegan to use the 
= facilities and their generous cooperation this work would not have 


been accomplished. The author is greatly indebted to Dr. Keulegan for the - 


interest and suggestions that he offered during the course of the experiment ce % 
‘uu to Messrs. L. Carpenter and W. Plummer for their assistance during the os 


setting up and running of the experiment. Finally, the author wishes to ex- 
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_ METHOD OF COMPUTATION ‘FOR STRUCTURAL DYNAMICS 


‘Nathan M. ‘Newmark lr, 


“This p paper is one the group of papers which were the bases for oral 


“The « entire group, , including inde se published in CIVIL ENGINEERING, will 


‘This paper describes a procedure of problems 
structural dynamics. . The method is capable of application to structures ll E 


placement ranging from linear elastic behavior through v various of 
¥ inelastic behavior or plastic response, up to failure. _ Any type of dynamic my 

_ loading such as that due to shock or impact, vibration, earthquake motion, or 

blast from a a nuclear weapon, c: can be considered. 

A method of numerical integration is described which for dedinae cases 
= for a relatively small number of i degrees of freedom is ‘suitable for use 

with desk calculators. However the ‘method is developed particularly for use 
‘with high- speed digital computers. Consideration is given to various types o of 
damping, and to nonlinear behavior. _A description is given of a method of 

treatment of elasto-plastic members in flexure, including the development of an 
yield hinges. By suitable means of application of the loading, and with the _ 

__ introduction « of enough damping to prevent indefinite oscillatory motion, the 

_ procedure can be used to determine the “static” behavior of a structure as it 


| through various degrees of inelastic behavior up to collapse. — Bi | 


anna for handling the problem of a structure having time-dependent a 
_ boundary boundary conditions, v which ch permits t the direct calculation of displacements a 


Note: Discussion open until December 1959, "Separate. payer 

_ mitted for the individual papers in this symposium, To extend the closing date on 

_ month, a written request must be filed with the Executive Secretary, ASCE. a 
_ 2094 is part of the copyrighted Journal of the Engineering Mechanics Division, —_ 

Proceedings of the American of Civil Vol. 85, No. EM 3, duly, 


} 
. 
Ap 
— 
. 

4 


ay The paper discusses the problem of structures with curved members and 

- describes in detail a method of dealing with arched or curved structures in 

which the members are not permitted to change in length. The treatment of 

‘such inextensible structures presents some difficulties in the general 
_ treatment but methods are available for modifying the general procedure by ; 

: _ introducing a series of constraint relations which reduces the number of de- 

grees of freedom permitted in the motion of the structure. 

ys The methods described herein have been used for the computation of the” 

; dynamic response of structures of various degrees of complexity including © 

is arches, domes, stiffened rings, framed structures, ¢ and simple spring- mass 

_ systems, subjected to various types of loading including nu nuclear bomb blast, 
—— easthagate foundation motions, random shock disturbances, wave action, and 
na impact ar and dynamic effects from moving vehicles. However, examples of 

‘ Ds _ The work on which this paper is based was supported at the University of 
-Dlinois by the Mechanics Branch of the Office of Naval 


general principles a and methods for are 
in this paper. The basic method of analysis is a general step- by- -step method 
of integration of the equations of motion, and is applicable to any ‘structure | § 
consisting of a group or series of concentrated or “lumped” masses supported 

oe ona deformable structure. re For convenience, the structure may | be considered 
to be a framework with joints or “nodes” at which forces may be applied, or — 
at which masses may be placed. _ Any finite number of degrees of f freedom may 

be considered, but it is essential in the procedure that the forces required to 

_ produce a pattern of deflection of the framework must be determinable when 

the deflections of the nodes at at which masses are are e specified od. It 

Because in an actual structure the mass is not enmay divided into separate 


eceasbbiaah parts connected d together by flexible elements, the structure which 
is analyzed is only an approximation to the actual structure. . However, a ee “a 
: i reasonable approximation can usually be made. For the actual structure, with 
F distributed masses, the number of degrees of freedom of motion of —o 
structure is infinitely large. The replacement of the distributed mass by a : 
ue number of concentrated masses reduces the number of degrees of freedom ‘ 

_ and affects the response of the structure in the higher modes. = “a 
The simplification of a structure for the purpose of making a dynamic 
me. analysis involves engineering judgment. One must select the essential ele- | 
ments of resistance of the structure so as to arrive at an adequate and reason- 

ably accurate expression of the resistance of the structure to motion. Also 

“ one must evaluate the masses of the structure and “lump” them together at 4 

_ places where they can be considered to act. In general the number of masses 


5 which h should be considered depends on the accuracy which one wishes “Tk 


; _ The structure may be made up of individual members connected together i 
joints would then be considered a as the nodes of the structure. 
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STRUCTURAL DYNAMICS 
However, the sinditeniiniiad an a one, either a one body or an as- 
semblage of plates and other elements, in which the nodes may only define m2 - 
points on the structure for convenience in the placement of loads and masses. i 
It is necessary that the nodes be so chosen that the resistance of the structure a- an 
to displacement of the nodes can be determinable, and that, if desired, the 
influence of inelastic behavior or plastic behavior can be taken into account. 
The structure is considered to be | supported at its foundations at nodes also. 
iy tes _ Although the method of analysis can take into account a situation in which Ye 
‘the masses change with time or with displacement, the procedure is presented 
herein for masses which remain constant. No major change in concept or in 
procedure is required to deal with the case of a variable mass-time relation- — 


ae Consider the plane structure shown in Fig. 1 which is made up of weight- : 


less but deformable elements supporting lumped masses. ‘The deformable ele- 
ments of the structure in this particular instance are shown as either beams Bs 
or bars which act under axial loading. However, much more complex structur-_ 

al types can be considered. . All the elements shown are deformable and conse- : wa 


quently each mass such as M, C, D, E, can move ove in both the vertical and hori- 


zontal direction. The points of support, such as at A and B, may also move | _ 


‘both vertically | and horizontally. ih 1 particular, these points may move in ——- = 


| have acting on it a force in any direction, or component f forces in the hori- | 
| The mass at M in Fig. 1 is shown removed trout the structure in Fig. 
Fig. 2a the point of attachment of the mass is indicated, and the positive di- 
rections of the resisting force exerted by the structure, R, and the dis- BEd. . 3 x 


placement of the structure, x, are indicated. In — 2b the mass is shown oe 2m 


have a vertical component of force acting on it also in which case it would | 3 4 
have a component of motion and a component of resistance in the vertical di- a 
rection, The sign convention that is chosen is determined by the arbitrary ¢t ; 
choice of the positive direction of the force P. The positive directions of x 

and of the resisting force R acting ¢ on the structure are taken a as the same as 


for P, and the | positive 


placement x are all the same. In general we will have for each pre i 
ponent of direction of motion of the mass a displacement, velocity, acceler- 
ation, resisting force, and applied force. The resisting forces R at any atin 
“of time are defined in such a way as to be that system of forces which are re- 4 
quired to pull the weightless deformable structure into a deflection configu- 
-Tation defined values of the x x at the same 
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Su "The ‘sign convention and notation a are > chosen so as to make it apparent that 


in modified form as R. If the forces P are applied very slowly there is only a 

a small acceleration and R is approximately equal to P. If the forces are ap- — 

plied quickly the difference between R and P can be very large. All the —_— 

stresses in the structure are Gcfined by the system of forces R. The xk 

structure can be analyzed statically for these forces. In general the forces R 

may continue to exist even after the forces P have dropped to zero. Similar- 

ly, R may be defined by the foundation displacements as well as by the de- ii 

flections of the masses, even when no force P acts on the masses. — aeal ie P 

. In general all of the ‘quantities so far defined are functions: both of position — 

on the structure and of time. If we consider a time t,, and designate the values 

of all of the various parameters including: displacements and forces at that 

S time with a subscript n, such as Ryn, our problem becomes that of defining the 
displacements Xn+1 and forces as well as the velocities and 

celerations @n41 at a time tn41, which differs from t, by the time interval h. ete 

In the argament which follows we shall deal with a single mass ina single — 

direction, but we might just as well deal with the whole set of masses =_ 7 : 

their possible displacements and designate the various directions with ad- eas i 

ditional m such as in However in the discussion which 
follows we shai. _rop the second subscript for convenience, remembering that a 
for each of the degrees of freedom | for each mass we have a a set of equations — W we, 
similar to the general set presented here. The derivation which immediately = 
follows is described in terms of a situation where damping does not exist, for =e 
convenience. . Later the procedure will be e revised to include damping forces. 

th general at any time, , (consequently or omit the subscript n), the Diiiitien 


4 It is that at time we know of the displacement. 
the velocity as well as the acceleration, but we know nothing about the situ- “ha es 


ation at time tp4;. Although there are . methods of numerical integration oahech 
permit us to make estimates, at least for small time intervals, of the dis- ye 


ing force during the time interval are not as sememeens as te mee described _ 

method presented here was derived by the ont first _ 

in Ref. 1. _ The relations which follow are given in terms of the acceleration 

at the end of the interval, an+l, although this is not in general known. A dis- 


earlier’ time, , these methods which do not take account of the change. in howell 


subsequently. Two parameters, Y and §, are introduced to indicate how much ie 

of the acceleration at the end of the interval enters into the relations for ve- 

locity and displacement at the end of the interval. The relations which are o 

(I + + 


| 
— 
” 
> 
ral 
— 
Is a spurious dam j 
ping introduced, proportional to the quantity 


It can seen if is s taken as a will 


involve a self-excited vibration arising solely from the numerical procedure, 
if Y is greater than 1/2 a positive damping is introduced which will 
reduce the magnitude of the response even without real damping in the © usa 
problem. Taking: account of tl the fact that Y must t equal WA we can rewrite 


e 
ban he 
In general unless B is 0 we may sroceet with our calculation as follows: 


wo Assume > values of the acceleration of each mass at the end of the inter- 
7 Compute the mass at the end of 
_ the interval from Eqs. (4) and (3), respectively. (Unless damping is 
.- present it is not necessary to compute the velocity at the end of the © 
interval until step (5) is completed.) 
oe a, (3) For the computed displacements at the end of the interval compute the 
_ resisting forces R which are required to hold the structural framework 
a (4) From Eq. (1) and the applied loads and resisting forces at the end of 
aaa the interval . recompute the acceleration at the end « of the interval. pe 
4 (5) Compare the derived acceleration with the assumed acceleration at the 
end of the interval. If these are the same the calculation is completed. § 
If these are different, repeat the calculation with a different value of as 
sumed acceleration. _ It will usually be best to use the derived value as 


the new acceleration for the end of the interval. rea | _ 


— The rate of convergence of the process toward equality of the derived and 
assumed accelerations is a function of the time interval h. Forasingle- _ 
degree-of- ~freedom system having a circular frequency of vibration 


BESTE 


ratio of the error in derived acceleration to the error in assumed acceler- 
ation, (where the error is the difference between a value and the “correct w 4 


=p=-8 : 6) 


ror in assumed acceleration 


s notation Eq. (5) can be rewritten as: 
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Now for convergence in a sequence of calculations the quantity p must b be 


less than 1. The critical value, for convergence, of the | time 
interval he can then be —— from Eq. (9) by 


Ly 


a Critical values of the convergence limit as a function of B are contained in a 


_ insure rapid enough convergence. If p = -0.32, the errors will be Teduced to 
x 7 e per cent of their original \ value in four steps or four r cycles of — 


_ ations affect the choice. The most important | of these is the matter of s aan 
For a complex system, it can be shown that the rate of convergence is _— - 
pendent upon the frequency or the period of the highest mode of the system. _ | 
: Consequently, the time interval used must be related to the shortest period of ed 
; _vibration, or the period in the highest mode of vibration, for the lumped <<a 
greater the number of masses into which a a system is broken down the 
_ shorter will be the permissible time int interval for numerical calculation of the 


Stability and Errors in Numerical Computation 


a ‘In order to study the stability of the numerical integration procedure, let et 
us consider the special case of a simple system, a mass with one-degree-of- we 
"freedom without external force acting on it. For such a condition, and for 
_ some initial displacement and velocity, the motion of the system should be a ao: ae ie 


with a circular = as given by the re- 


which a is the spring constant and M the mass. The: tri between 


he acaeration and the displacement is g 


use of th the symbol as defined 


a the lower values of 8 would be best to use. However other consider- i | 


; we can derive a ae chen relating 1 the values of three successive 
of the equation in in general terms is: 


the relations between finite and derivatives, 
is = be seen that the last term on the left of this equation corresponds to a 
factor times the velocity of the system, and consequently can be iateepeeted 
& a viscous damping term even though the system was defined as having | 
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— 
no comping. ‘This sp spurious last term can eliminated by the choice 


‘The g general solution of the finite difference equation , Eq. (15), can lis 
written in the case where the quantity a2 is less than 4. In this case, define a 


17 

This c: can be of a pseudo period and an initial dis- 
placement x, and a B which is) the same for ‘form as velocity. ‘The 
itt 


: This may be | compared | with the exact solution, x, which is given ren by Eq. (19): 


cos 2nt/T + sin 


It ca can be seen that the approximate solution, Eq. (18), is to the 

exact solution and gives precisely the same response for an initial ge cid 
a placement, but gives a different period from that of the actual —om The 


value of the pseudo period Ts ; 


ine 
“The fol may be used for a simpler 4 definition of the 
relative value of the pseudo and the real period of vibration: ea 


8 2724 - 1208 + 7208 


Ther response of the sysiem to an initial velocity is given by the nea 
4 er in Eq. (18). The reiationship between this response and the true re- 
: cs sponse to an initial velocity, as shown n by the second term in Eq. (1! 9), is indi- 


Bi is exactly to 1/4, the | velocity response is co correct but 


‘it t is different from 1/4 there is an incorrect maximum velocity response. 
a = Values of the errors in the period and of the errors in maximum response 


and for a range in values of h/T. There is also given in Table 2c the rate | of 


convergence for the corresponding tabular entries. Fora system with 
‘number of degrees of freedom, the limits are expressed in terms of the 


— 


= es When «2 > 4, the solution of Eq. (15) oscillates without bounds, and the 


to an initial velocity are given in Tables 2a and b for several valuesof 8 


calculation does not yield results even in with the | exact 
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CONVERGENCE AND STABILITY. LIMITS 


Convergence Limit 0.450 0. 389 0. 
Stability Limit, h/T 0.389 0-450 0.551 in 


pressed | of the time interval relation between a 6 in 


er 


re 


which can be simplified to tl form: 


Bees Table 1, it can be seen that for values eh greater than 7 if the 
time interval is chosen for convergence the numerical procedure will always ; 
Fe stable. However, for values of 8 less than 1/8, convergence does not in- 
sure stability. Lack of stability gives no warning of difficulty, but introduces a 
Ba spurious increasing oscillation into a system which may be in oscillation — 


P anyway. Therefore an inexperienced computer may not recognize the diffi- 


parent to an experienced computer, - Consequently it appears that unless other 
steps are taken to insure stability, one should limit the time ngewe all e 


A method very much similar to that described here for § = 0 has been dis- ra 
cussed in Ref. 2. A method corresponding in many respects to that 
_B=1/12 has been given in Ref. 3. _ However, the general treatment previously — 
presented is different from that given here, particularly in the treatment ot 
the starting of the motion. A method similar t to that for 8 = 1/4 was first — ; 
presented by S. Timoshenko in Ref. 4. . However, he did not carry the pro- 
cedure to the point of generalizing it for other than simple yn uadinallls 
freedom nor or did he develop the ® conditions stabilityand 
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EFFECTS OF LENGTH OF INTERVAL ON ERRORS DUE T DUE TO NUMERICAL PROCEDURE 


Relative Errore tn Perko 
076 -0.006 0. 028 059 0.12) 
. 70-015, 0.038 087 oO. 


(b) Relative Errors” “im Maximum: Response to an Initial 


0.025 0.017 
0.20900 
0.614 


inf. 


0. 0.008 (0.012 0. 
0. 033 0.066 


500 a 750° 


corr 


Pps to a step function with a uniform value hyve to the initial value for th 


first half of the time interval and a uniform value equal to the final value for | 


the second half of the time 
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; 4 ail valk poe four values of it is possible to define 
possible for vaiues, al. Three of these are 
f acceleration in the time interval. 
, | if - in Fig. 3. It appears that a choice of 8 = 1/6 cmmmespo | 
n in Fig. 3. It appears that a cho 
shown in Fig. 3. 4 
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< the beginning and end of the time interval with each double pulse consisting 
of a part equal to 1/2 of the acceleration times the time interval, one occur- 
= just before the end of the — isterval and the other just after the 
beginning of the next interval. 

Treatment bars ‘ 

Dampin ing 

forces may be proportional to nny or to some power of the velocity, or 


a lute”, and depend on the particular motion or velocity of the mass with re- 
spect to ground, or “relative” and depend on the motion or velocity of the 
with respect to some other point on the structure. A structure with several 

: ‘points having damping forces is shown in Fig. 4. If masses M, and Mg are — 
different from ‘zero, then the damping forces can be applied to the masses, as 
indicated on the 1 right- -hand side of the figure, and the calculation is made in : 
After the acceleration at the end of the interval is assumed, the velocity at 

_ the end of the interval and the displacement are computed from Eqs. (4) and 
(3). The complete motion of the system is now defined at the end of the inter- 

val, and regardless of the kind of damping, the damping force is determinable, 

- With the | damping force placed on the mass in the proper direction, the ac- 

_ celeration at the end of the interval is now computed taking into account the 

4 ‘damping force as well as the applied force and the resisting force. . If this is: 
Giflerent from the assumed acceleration the calculation is repeated. Sie 


is 
\ 


| 


) structure are computed from the , displacement of the the ‘structure at the end of 
= interval. ~The damping forces acting on the masses are computed from 
the displacements or velocities of the masses in accordance with the > damping 


a It can be seen that in the case of denabia ai value of B = 0 | presents no real 
_ advantages over any other choice of 8 because of the fact that the accelerati® 


oO 


locity o1 or other r parameters ; which . determine the damping forces at the end of 

i For hone ¢ cases where there are no masses at points w where damping force 

_ or external forces are applied, the situation is more complex. “Fig. 4 illus-— 

- trates this by applying forces to the two zero masses, My and Mg. It is necet 

sary in these cases that the resisting force be equal to the algebraic sum of 

the external load and the damping force, and this resisting force is then ap- 

_ plied to the structure and except for its dependence on the damping force re- 
mains at this particular value, regardless of the displacements. However, it 
La is now necessary to compute the accelerations and velocities of the point on 


aM the structure in order to find the damping force. This can be done in the a 


Get 
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following way for the points where the massis zero: = = 


Assume a value of acceleration of the mass 3s point even the mass 


Compute the ant | displacement | the point and d determine the 
damping 
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} modal behavior for structures which are not elastic, it is not desirable nor 


STRUCTURAL Dy DYNAMICS 


(3) Now apply the net resistance, corresponding | 7" the difference between uw 


ou the applied force and the damping force, to the structure and determine ~ 
the displacements of the structure at ‘all points when the qremeSes dis- 
_ placements are put in at the points where masses exist. = ~ 
_ (4) Compute the acceleration at the end of the interval that is required to 
give the displacement determined in the preceding step. 
(5) Compare the acceleration so sO settee with the one initially assumed 


‘This s procedure is is considerably m ‘complex than that which is used 


“masses exist at all nodes. Consequently, it may be | better to put in an actual — 
: mall mass than to make the mass zero. _Sowwer, if damping forces are not 
relatively simple to 
handle the problem d directly. 


Recommended Time Interval and Choice of 8 
su 


In any set of calculations in which an error in any one step roe all = 
bsequent steps incorrect, it is desirable to have a self checking procedure. 
- Consequently the procedure described herein works best where it is used with 
= aconverging system of approximations because under these conditions the 
| calculations in any one time interval are repeated several times with slightly 
. different values of the numbers involved. A close agreement between the re- _ 
_ sults of successive calculations is in general a ‘sufficient check on the accu- 
, racy of the numerical work when the work is done on a desk calculator. (Suc * 
rs a check may not be necessary on a high- speed computer.) It canbe seen, _ 
: therefore, that a value of B = B = 0, , where there is no damping, is not necessarily 
-a good choice because a separate check will be required in these casestoin- 


ae Studies of the effect of damping and of negative spring constants suchas | 
- those corresponding to a decrease in load with increase in displacement, indi- 
‘ cate that better results are obtained with values of B in the range from 1/6 a 
_ to 1/4 than in the range below 1/8. In general, with a time interval of the _ 
_ order of 1/5 to 1/6 of the shortest natural period of vibration, the rate of con- 
is _ vergence will be rapid enough for all practical purposes, and the errors will . ¥ 
be small enough to be tolerable for every combination of damping or negative 


_As a structure goes into the inelastic range, in general the periods of vi- o 7 


_ bration all become longer and the shortest period becomes longer as well. Be oes 
af For purely plastic resistance at the limit condition of an elasto-plastic _ ed 
_ structure, the period is infinitely long. Consequently the time interval can be- ; 
generally desirable to increase the time interval in accordance with the change 
in the structural rigidity as the structure becomes inelastic. However, it is 
Rot always | convenient to compute the shortest period of vibration for oo. wi 
Bey! which goes partly into the plastic range. Consequently it is desir- ; 
able to establish the conditions which govern the choice of time interval on Md 
/ Some measurable behavior of the structure which is a natural function of the a 
_ method of calculation. Since in general it is not practicable to consider __ ets 


J come considerably longer as plastic action develops in the structure. It is ae . 


—— to separate the individual ‘modes of action of the structure. ‘gue 


the 


an the use of a value of B greater than 1/8, the more _——e of the ca 


| 
| 
re- 
SE 
oc! 
minabl 
this is 
end of 
from | 
1amping 
no real = 
ever, 
ne the 


convergence will be an an adequate criterion for the time interval. (2) 
_ ly, one can establish a rate of convergence based on the number of iterations — q 
which ‘it is desirable to make in a time | interval and then examine | e the rate of 7 


ample, if it is desired to nae three significant figure accuracy in dis- - ie F 
_ Placement velocity, three significant figures i in the acceleration are also 
desirable. A rate of convergence such that the error is less than one part in 


‘time interval to natural period of the order of 0.10 for B = = 1/4 will nth this 
— rate of convergence. Similarly, a time interval slightly longer will be ade- 
quate for 8 = 1/6. ‘However, it appears that these time intervals will introduce the di 
_ errors of the order of 2 per cent to 3 per cent in the period. If this is not ads 
missible, then a faster rate of convergence to the 
ss Let us say for example that a time interval consistent with a rate of c con- 
a ergence of 10 per cent is desirable and that three cycles of iteration will 
_ normally be considered convenient so as to ) bring the results” to an error 4 
7 less than one part in a thousand at the end of the third cycle. ‘Then a fourth 
_ cycle will verify that this is in fact the case. Under these conditions, then, — 


_ We can assume a time interval and run through t the calculations | several times. 


too long a time interval and we can shorten the length. Since the rate of con- 
_ vergence is generally a function of the square ¢ of the time interval, the next . 
- estimate can usually be decided on from the course of the preceding x calcu- 
lations. If the time interval is made smaller in accordance with the estimate — 
, ia of convergence, , and if convergence is now in fact obtained with about th three to 
“ag 7 four cycles of iteration, then this time interval can be used in subsequent " 


steps. If however the time interval chosen leads to convergence in only two. 

7, = or three cycles, we have probably taken too short a time interval and in the © 

i /* cycle in the calculations we may choose to use a longer time interval. ae ent: 
_ With this type of procedure, and with continual reexamination of the time | 

+ x5 _ interval in terms of the rate of convergence, one can take account of the 

E 4 change in the characteristics of the structure without loss in accuracy and 


ma 
supplementary calculations. The procedure described above can of fr 


readily be programmed for a high-speed digital computer. However, an upper 
= on the time interval is usually desirable to avoid difficulty in terms of Pa 


a 

, and with deformation in all of the members, it is a relatively by c 
simple matter to perform the calculation for the resisting forces if the dis- stru 

_ placements ar are given. . All that is required is to take the set of displacements ® plac 
. at a given time instant, pull the joints of the structure into the corresponding © foll 
er compute the > axial forces in the | members from the changes in in a 7 
aa length between the nodes, and compute the moments in the members by Be _ 


rocess which involves two stages as follows: 


Consider the joints locked against and determine the fixed end 
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E (2) Relax ses joints by permitting them to rotate. Either moment distri- " 
bution or some other technique of calculation can be used to determine — 

‘The’ shears and axial forces can now be computed in all the and 
‘one can compute the horizontal and vertical components of the forces on the 
pins ' which hold the structure in the displaced configuration. . These forces, = 
reversed in direction, are the resisting forces. acting on the structure. 
_ If the structure is not elastic, it is possible to make the calculation for < any 
known load- deformation or axial-displacement deformation or curvature- 9 
moment relationship. . The technique of the calculations is unchanged. How- 

ever, for nonlinear behavior it is essential that one take account properly of oa 
the direction of loading in a member or element so as to be sure that the  . 
ment is continuing to deform in the same direction with the appropriate re- 
duced stiffness, and to determine when it is recovering or unloading with the ay 
appropriate increased stiffness. _A test of the direction of relative deformation 
can readily be made to determine this. 
_ Ina structure such as shown in Fig. 5, where masses might be considered r 
placed at all the joints of the structure, it is usually neither convenient nor _ a 
desirable to consider the most general situation in which the masses can be p 
displaced horizontally as well as vertically. The reason for this is that the , 
period of vibration of f the masses in the horizontal direction is very short 4 
compared with the period of vibration in the vertical direction. If one were wee 
= to solve the problem with the general procedure, one would have to use anex- 


tremely short time interval which would make the calculation tedious wa 


th such cases it is usually desirable to consider a restricted type of defor-_ 

® mation of the structure to reduce the number of degrees of freedom. This — a 

® can be done by connecting the masses to the structure by means of vertical ; 

3 links as shown in Fig. 5, and specifying that only vertical forces can act on the 

= masses which in turn are considered to move only vertically. The links an a 


considered to remain vertical by making them long enough so that the hori- _ 7 
zontal motions of the truss joints do not introduce an angle into the link - 
tween the trussjointandthe mass. 

| Ih the structure shown in Fig. 5 the number of degrees of freedom with a 
mass at each node would be 20. However for the system shown only 5 degrees 
of freedom are needed. Both the difficulty in the calculation and the number | 
of repetitions of calculations for a given duration « of motion will be greatly in- — 
creased if the 20 degree-of-freedom structure is used. - However, some diffi- 
culty in computing the resistance functions in Fig. 5 is encountered because : 3 
it is not possible directly to determine t the changes in length of all the ie 4 ~ > 
members from o only the vertical component of apacement of the lower r chord 7 


In cases of this sort, it is possible to proceed in a . slightly different fashion ges 
by determining either directly, or by inversion of the influence matrix for oad fe 
structure, the set of forces required to produce the given set of vertical ai 
placements. In the general case, , one can summarize the calculations as 


br @ Compute the vertical displacements of each of | the masses for individual 


values obtained is as the influence for deflection 
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the forces acting « on each of the masses for a unit roots Aeneid of ll 

~ of the masses in turn. If this inversion is performed once, the matrix _ 

ae values can be applied to any set of displacements to determine -_ re- 

7 (3) The coefficients in the influence m matrix will change when any elements py. 

in the structure become plastic. Consequently the calculation of the 


> 


influence matrix and the inversion of the matrix to obtain the stiffness > 
— iF matrix must be performed for each change in stiffness that occurs dur- 2 
— oe a... the history of the structure. It may be desirable to make the calcu- 
a lations by a “relaxation” procedure for more complex structures of this 
a 3 _ However, in principle the structure can always be analyzed to find 
the for 


it might be pointed out that for a a beam of the same span as in Fig. 5, wi ith 
Ba series of masses acting on it, it is not hecessary to use the influence matrix 
Sand invert it to find the stiffness matrix. The structure can be gn di- 


ihese deformations to find the shears in the members and the consequent — ei 
forces acting on the masses. ~ Although in such a case there would be theoreti- 
Pw cally 10 degrees of freedom, 5 of the degrees of freedom are associated with 
j a ‘@horizontal axial deformations in the beam, and these are not “coupled” with 
> the vertical deformations because of the nature of the system. This obser- ae 
| _ §§vation is true only for a straight beam. In a later ERE the procedure Beat 
_ The structure in Fig. 5 is one of several apecial cases s which require a eC 
slightly modified form of treatment. Another case concerns any structure R= 
with constraints on the deformation of some of the members. Such a structure — : 
might be a polygonal arch which is considered to deform only in flexure and = q 


which is not considered to have any axial deformation of the arched members. 
In this case also the number of degrees o of freedom are reduced, but the Bsns 
reatment must be somewhat different from that for the truss because 
al umber of constraints m2 may reduce thi the e number of degrees of freedom even 4 
— a ow the number corresponding to a mass at each node prescribed to move in > ‘ 
7 4] 2 particular direction. Further discussion of this topic is contained in the _ Ao 


A general procedure for such structures can be discussed. In general, , one 
ecan handle the problem in the following way when damping forces are not _ i 


present. . The modification for damping is readily made as has been | 

Assume values of acceleration only at the joints where masses exist 
and in the direction in which the masses 

be —@) Compute the displacements of the masses, but not the ye displacements rot ~ 

« Let the joints without masses be free to deflect but meres to a joints 


xternal 


Both Figs. 4 and 5 show wi with zero masses at some of the nodes. 


in 


From these influence coefficients, compute the set of forces on the 

a structure for given values of the displacements. The calculationin- 

volves solving a set of simultaneous equations for each set of valuesof 

This can be svstematized by “inverting” the influence _ i 

4 
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y Elasto-Plastic F 


_ forces are zero, no ote and no loads are applied at the massles; 


(4) We now have a a system with prescribed deflections at some points only 
a. and with prescribed loads at other points. It may be pointed out that at 
are some particular joint for example there may be a prescribed vertical 

_ displacement anda prescribed zero horizontal load, if we consider that 

id — vertical displacement of the mass is to be considered at that poin 


the pronto of motion. In general this can be done by : an influence . 
—— gece in which we apply unit forces at the masses whi 


“ _ which in n effect says that the product of each unknown reaction force at 
ri each mass multiplied by the influence coefficient for deflection at each 
ie of the masses, and summed for each of the masses individually, leads 
be to a deflection at each mass equal to the prescribed deflection minus > 
the deflection in that direction at that mass due to the prescribed force 

acting at the nodes where no masses exist. 
 o The solution of this set of equations yields the desired results. We ma 
_ now have to complete the calculation by determining the displacements 
: ibe at the massless nodes of the structure if these are required for any — 
eos af other purpose. In most instances they will not be needed and need not 
_bé determined. Where damping forces exist, however, we will need 
these displacements so that we can determine the velocities and accele 


forces on these quantities. 


fi Non-Elastic 
General Comments 
| 
_ For trusses the problem of inelastic behavior of the members isa relativg 


Pi ly simple and straightforward one. One n need take into account only the cha ing 
7 oe in deformation of the members in each time interval to determine whether hee 


general « concept for framed members is the s same, there are some difficulties: 


= resisting in the structure are directly Although 


The method for « dealing with elasto-plastic behavior in frames is describe 
_ here . Consider a general framework consisting of members acting | in flexuré 
as shown in Fig. 6a. At some stage in the deformation, the moment at the | 

~ ends of the members shown a as black dots in ' the ' figure h has reached the plasti 

iad limit moment. One or several of the members meeting at a joint may be lo val 

sen . a general enough to consider members which are not loaded internally 


in which the moments a arise from a set of of th 


action is or the member is unloading 
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sslesy an interior load applied dynamically, Then we will in effect | have a prescribed 
displacement of joint E which can be considered to be an external hinged sup- 4 
s only port as shown. Consequently, our method of treatment is general eaenasiis deal 
that at with internal loads as well as loads applied at the joints. om 
‘tical eu At some time t, the structure has been analyzed with the result that we 

er that have determined the plastic hinges which exist at that time and, as a ae 
t point quence, , we have determined the velocities, _ displacements, and resisting forces 
sses in) in the structure. We can also determine from the final configuration and Caan 
nce & moments in the structure at that time the rotations of the joints 6; and the ro- 
s whi tations at the plastic hinges . if we have analyzed the structure by moment 


hese § distribution, the fixed end moments can be computed and distributed and from — 
uations a comparison of the fixed end moments and final moments the rotations at g ie 
rce at each end of each member can be determined. . The rotation of ajointisthe | 


ut each same as the rotation of the ends of all members meeting at that joint which do 
leads § not have plastic hinges. The “kink” or relative angle between a member hav- 
ninus | ing a plastic hinge at the te and the angle of the joint is measured by the : 
d force difference between 6, and 6;. The kink angle is of importance in sonaiiadlll 
| whether the member is or unloading at that joint. 

considerably simpler to deal with the structure in each interval 
ements _ making calculations of the change in the structure. In other words, one works — 
any — with the change in displacements during the interval as defining the fixed end es 
ed not — moments at the end of the interval and the final moments so computed are the 
need _ changes in moment to be added to the original moments in the various nee 
| accele 


members. In this process one can deal with the elasto-plastic hinges as actu- 


ping _al hinges if the structure is loading during the interval and the hinges will = 
al removed if the structure is unloading during the interval. 
If the total moment, including the plus the increment 
iol i moment during the interval, at the end of any member exceeds the 
hinge moment, then a hinge must be placed in that member. 
relat ' @ ) If a hinge exists in a member and the increment in kink angle in the 

a al Cae member at the joint at the end of the interval is in the same | direction | 
as the preceding total kink, then the hinge remains. 
laving | (3) If the kink angle in what was originally a plastic hinge should decrease, 
orces Of) The structure is modified by adding, keeping, or removing the ng 


¢" ee _ Only in unusual cases can there be plastic hinges in all members meeting 


‘a 7 - ata joint. Such a case is shown at D in Fig. 6. In such an instance the joint — 

rotation can be assumed to remain at its value at the time the last members ies 
at the joint become plastic. (It is clear that in such a joint at least two i 
members must become plastic finally, and at least two members must simul- ei 


lescribe 
n flexurt 


eg - member F in which both ends develop plastic hinges. Such a member can be 
scan considered to have end rotations the same as those which existed prior to the — 
nts of the development of the last hinge in the member. 
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na structure as the loading on the structure is applied, moreorless = ~~ ~~ SS 


“statically”, r applied ‘released in various ways. ‘This ca can 1 be done very 
a simply by specifying a relatively slow rate of application of the loading, corre- 


vibration in the fundamental mode. ‘usually « convenient to introduce 
- damping so as to avoid oscillatory motion during or following the ening 


are not concerned with the dynamic res response, § so as to make the - fundamental 4 


_ — if the  iinidiaioen or supports of a structure move with time no 
~~ change is required in the method of calculation and the general procedures de- 
scribed herein are directly applicable. However, there are instances where 
it is desirable, , or convenient, to deal with the motion of the elements of ah 
structure relative to the base or foundation of the structure, rather than in 3 
absolute terms. Such situations arise commonly in earthquake motions. _ The 
~ method described herein is useful primarily in those cases where the “wid 
- dation moves as a single unit and where motion of the foundation of the any 
_ structure with all masses 5 equal to zero would not introduce stresses in the 
_ Consider the structure shown in Fig. 7 where the base can move as a unit 
E in the horizontal direction. Let the motion of the base be defined asa function 
“of time by the quantity y- ‘The ‘motions of any of the masses in the structure ; 
are defined by x. Let us assume that the axial deformations of the members| 
can be neglected and that we have only flexural deformations to consider. In 


~ such cases, we can lump the masses at each floor level at one point as shown. § 
more general case offers no difficulties, however. & 
a a atl If now we —" to the structure with masses on it a force F, in ete a 


id din the structure and the foundation will move as a unit aids no relative dis- ¢ 


ex an now we apply to the structure c considered supported o on a fixed foundation . 


force Pp,  datined by the 


ah 


a time-dependent boundary motion into one in which we have external forces 


we have the from one which which there are are external forces and 


and no boundary motion. The for th the | Pp will 


ay 


it can soca, of course, that th the total motions of the ons hei 


i of the relative motions u plus the foundation motions y, and the total loads ap-) 


plied are the external forces P. . Consequently, the procedure is valid. 


a For more complex systems, one must modify this procedure by taking into 


; ‘te account the displacements at each of the several points of support, and defin- 
; ing the quantity y for the motion of the masses as being the motion at that as 
- = a ae structure consistent with the foundation motions. In the 
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case se ofa a ‘for ex exam] mo one and hinged at 
_ other where the hinged end is constrained to move vertically in accordance > 
with some time-dependent condition, the displacements o! of the masses along = 
the length of the beam are proportional to the end displacement and are equal 
to that end displacement multiplied by the ratio of the distance from the re 

_ Support 1 to th the 1 mass divided by the length of the beam. _ The same method of 


“of the beam, relative to boundaries fixed in position. 
In this case, and in other cases where no stresses are introduced inthe _ 
weightless beam by the boundary motions, the entire solution is — , 
forward and all of the stresses can be obtained from the modified system. 
However, in more complex cases where the boundary motions may introduce — 


‘stresses in . the massless system, the general principles are applicable but the f 


ay ries and with modified forces. In such instances, the problem of plastic be- 
havior introduces complexities which may be unwieldy and it is not desirable 
to work with the replacement dsp In all cases it is 8 possible to work with 


> Curved Arc 


"General Procedure 2. 


- Consider the two structures ‘shown in Fig. 8. In Fig. 8a, a beam is shown 


‘with two masses. These masses may move vertically with consequent flexure 


of the beam or horizontally with | consequent axial deformation. There is no 
_ coupling between these sets of motions. There are two degrees of freedom in — 
each set or type of motion and four degrees of freedom altogether. The verti- 
cal motions can be dealt with as if ina two- -degree-of-freedom system and the J 
axial motions can be dealt with also as ina two-degree-of- -freedom system 
with neither being considered to have influence on the other. However, in the 
structure shown in Fig. 8b, neither pure vertical nor pure horizontal motion ) 
of each mass is possible. In general, all four modes of motion or all four de- 
«grees of freedom are coupled, and the structure requires treatment if 
_ However, in the structure shown in Fig. 8b, if it is desired to avoid m 
deformation of the bars, only one mode of motion is possible, with the ahanel 
_to the left moving down to the right and the mass to the right moving up and to 
the left. The system has only one degree of freedom, and can be analyzed as 


Consider the structure ire shown in Fig. 9a. Here. are interior 
masses and if the members are considered inextensible there will be two de- 
ae grees of freedom. We ve arrive at this number by considering the total number 
degrees of freedom of the system if the members were capable of defor- 
mation, namely 6, and subtracting from this number the of 
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stresses must then be determined by adding the stresses in he massless 
"as _——~@ system with boundary motions to the stresses in the system with fixed bounda- pe 
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= in terms of the relative displacement of joint C compared with aad te 
joint B and the s slope of the member BC. For small deformations, where 5x i 
and dy represent the horizontal and the vertical motion of point C with respect 
to point B, and where 6V and 6H represent respectively the vertical differ-— 


ence in elevation between C and B and the horizontal distance between Cc and 


can be expressed as follows: | 
This equation can be iieiaaeckias as ici that the deformation is such that ; 


2 

‘point C moves perpendicularly to member BC. 
A similar constraint relation can be written for each of the bars. If the i 

total deformations are expressed in terms of the increments in the defor- se 

_ mation, one can find four relations among the three horizontal and three verti- 

displacements of point B, C, and D. These will be somewhat more com- 4 

¥. her plex but the relations can be written in any case. _ The number of independent pl 
; b relations which it is possibile to write in this case will be equal to the number 


of constraint conditions, namely four. - With four relations ‘among six x dis- 


Bk: maining. These can be chosen as any ‘two components of displacement, “a 
_ the calculations carried through in the usual way, but with only two masses 
nal 7 and two directions of motion considered. . By the use of the constraint relation- . 
Mi ships, when the two deflections which are independent are assumed the re- _ ‘te 
‘maining deflections can be readily computed, and the entire ee: of defor- 4 
_ However, it is possible to develop a more direct PRE ck for determining = 
_ the relationships for the resisting force components a and for the displacements. 
This can be done by noting that there are certain patterns s of loading on the > s 
‘structure which produce only axial stresses in the members. For these load- § : 
ing components no deformations of f the arch can take place. It is obvious that ¢ 
: one type of loading would correspond to axial forces directed along the length 4 
of any one of the members. This can be expressed also in terms of the verti- 
eal and horizontal components at the two ends of of the seated such that the 
lies along the length of the member. wie 
__ By considering four component loadings ancients to loadings produc- 
‘ing axial force in only one member at a time, one arrives at four * independent 
— load conditions which produce no. no displacement in the arch in Fig. 9. These a : 
_ can be combined in various ways. One obvious combination consists ina 
pattern of only vertical forces at the interior nodes. These of course would 
be resisted by vertical and horizontal reactions at the ends but these do not 
_ enter into the picture if the ends do not move. This pattern of loading is shown 
in the upper part of Fig. 10. A second pattern of loading in the lower part of 
v. ie ‘Fig. 10 | consists of the two sets of loadings which produce tension only in in the 
a) eater members of the arch. If the structure is symmetrical, and if the loads 
: ee B and D are made equal, then the loading pattern will be symmetrical. ‘The > 
—_— pattern in the . upper part of Fig. 10 will always be symmetrical for a 
symmetrical structure. These two loadings represent the two 
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symmetrical — producing no displacement. italia: one can arrive * 
at a loading consisting of a downward load at C and equal horizontal loads with — ‘ 
no upward load at joints B and D. There are two other independent loadings j 
_ which are anti-symmetrical. From combinations of the four loadings, various — 

| types of loading patterns can be derived. 
. ‘The advantage of having the constraint conditions stated in this way is is that 7 
} now we can add any of these types of loadings to any set of resistances or of Bs 
_ external loadings so as to preserve the number of degrees of freedom of the ~ - 
_ structure at the proper number, in this particular case two. In other ro 
2 after the two independent yg army are decided on such as the horizontal : 


_ make the accelerations in the various Seadiens of all of the 1 masses es con- —_ = &§ 

Consistency in the vertical displacements during the requires con- 

_ sistency in the vertical accelerations at the end of the interval. Since the dis- 

A placements at the beginning of the interval have been made consistent, and 
since the velocities and accelerations at the beginning of the interval are con- 
sistent, we need only insure that the accelerations at the end of the interval — 
are properly consistent. In order for these to be properly consistent, they 

- must obey the same relationship that the displacements obey. In other words, 

_ we have the same relationships among the accelerations in the various di- | 

‘ rections at the joints or nodes as we have constraint relations corresponding 
. to zero extension of the members. It is a simple matter to adjust the resist- 

_ ing forces by adding the proper components of loadings so as to make the = 

acceleration at all of the masses consistent. — 

by-product of the results shown in Fig. 10, one can apply the 

_ of virtual work to find relationships among the displacements at the various | te: .- 

_ joints. For any pattern of loading which produces no displacement, the = i 
_ product of the loadings in the pattern times the displacements in any set of 
: ‘conatatest displacements of the structure in the direction of the loadings in ba! oe 
- the pattern must be zero. This follows from the principle of virtual work ~ 
‘i rectly. - In the case of an arch the vertices of which lie ina parabola, asin ot 
5 Fig. 9, the loadings at the three vertices, as shown inthe upper partof = 
Fig. 10, will be equal. This implies that the sum of the vertical deflections at = : 
Ss C, and D, taken as positive downward, must be zero. A similar relation- al a i 
ship can be arrived at from the loadings in the lower part of Fig. 10 = from / 
each of the load conditions considered. 


5 


c analysis of an arch which is inextensible. Consider the arch shown in Fig. 11 ba 

_ which is the same as in Fig. 9 except that the masses are carried by ai ae 
_ links, which remain vertical, in the same manner as in the truss in Fig. 5. as 

: ‘Here n now we constrain the masses to move only in the vertical direction. a ia 
_ There are then only three degrees of freedom for the structure. If the BP yee 

: structure were extensible, we have reduced the number of degrees of freedom o 

" from 6 to 3, but we have lost the horizontal motion of the masses. If we know 

a iat these are unimportant or feel that they might t be, or if we vam to — a. 
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Analysis of a an “Arch” 
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- same fashion that we solved the problem in ‘Fig. 5 earlier. We need ni now y find 
only the influence matrix and the corresponding stiffness matrix in order to 
however, the arch is inextensible, there are ‘really only two degrees of 
freedom of the structure and we have a system which apparently has three. 
However, , by use of the fact that for the arch shown a set of loadings consisting 
fs of equal vertical loads at each of the three joints will produce no deformation, 
we arrive at the condition that the sum of the three vertical displacements - 
r _ must be zero. | Then we could use this as a constraint condition which will 
permit us to work with only two independent vertical Gapfincements. We can 
_ define the third vertical displacement in terms of the other two, and arrive at A 
a consistent set of forces by adding to any set that we compute at the two inde- 
pendent joints a proper combination of three equal vertical loads so as 
_ make the vertical accelerations at the three points consistent. Consistency — 
in the vertical accelerations means consistency in the increments in de-— “—— 
flections, which means that the vertical accelerations must also obey the rule 
-* that the sum of the vertical | accelerations a at points 1, aa. baal 3 as sho shown in emg 
It can be seen that in general for an arch with N nodes, the number of de- = - . 
grees of freedom n for an extensible arch will be 2N, and f for an inextensible — me &g 
N-1. The number degrees of freedom for vertical motion at the 


cee If the ends of the arch are free to move, _ additional degrees of freedom ae é 


the extensible or inextensible case ina simple manner. 
"ever, ‘more complexities are introduced in the case of a complete ring. Bite tte 
lia _ Consider, for example, the hexagonal ring in Fig. 12. There are six mass- 
es at the nodes of the ring. For other shapes or number of sides, the re- ee 
_lationships can readily be determined. In this particular case, it is clear that $ 
_ the number of degrees of freedom for an extensible ring is 12. However three a 
a ‘a these are the rigid body motions for an object in a plane and produce no in- 
) ternal stresses. In the general extensible case of deformation, with two de- 
3, grees of freedom for each mass, , one would compute ordinarily 12 dis- Bi 9 A, a 
placements, one in each of the two directions at each of the 6 masses. Rai | ' 
_ Statically required forces needed for consistent accelerations will automatical- 
ly be achieved, and the system can be treated as if it had in fact * | degrees of Aes 
In the case of the inextensible ring, there are in the case of the hexagon 
_ six constraint relations corresponding to the six bars. ‘This means that — 
- are essentially only 3 degrees of freedom for the system rather than 9 for the | 
inextensible case. _ ‘Here again the 3 rigid body motions correspond to ron, Hy 
_ stress in the system. However the system can be treated as if it had in fact Pica 
ie six degrees of freedom since the rigid body motions will automatically be ae 
achieved. The six constraint relations must be used to insure consistent _ ie 
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4 duly, 1 
Miscellaneous Problems 
Space Frames and Space Structures be 


‘Space frameworks or space structures of a more complex character than 


) 2 ples. Fora a structure in : space, there are in general 3 degrees of freedom tad 
a 4 each mass at a node, and consequently 3 components of external force and re- 
" sisting force which must be considered. | This introduces complication only a 
in the order of magnitude of the calculations which are involved. ca a ie 
For | in which axial loading produces secondary effects because 
of the deflections of the parts of the structure, thereby increasing the 
"moments and the deflections for a given loading, , the calculations are more a 
; _ complex. In general, the resisting forces for a given amount of prescribed 
: a deflection are reduced, in about the same ratio that the deflections for a given 
loading are increased when buckling loads are present. ‘This reduction in re- 
_ - sisting force means generally a greater acceleration and a consequent greater 


_ dynamic displacement in a structure in which buckling tendencies exist. _ These | { 


tendencies can be taken into account but the methods for doing so become rela- 
_ tively unwieldy. One can compute the resisting forces on the assumption that 
a buckling tendencies are negligible, and then for these resisting forces de- 
: ies termine the relative deflections corresponding to the secondary effects of the 
forces produced by the deflection of the structure. Then by reducing the re- 
; . sisting forces in a more or less arbitrary manner, one can arrive at a set of 
f resistances which would account for the prescribed deflections. a Although this 
_ method is not entirely satisfactory, it does give a means of taking into account 
approximately the buckling tendencies for complex structures. For simple 
_ structures, the problems can ‘sometimes be handled directly. The real diffi- : 
~ culty in the case of buckling problems is that buckling is a nonlinear effect _ 
= the methods of calculation used to handle the calculations in com interval . 
in the process depend on linearity, within that interval at least. — 
“4 
a . _N. M. Newmark, “Computation of Dynamic Structural Response in the a 
pet e Approaching Failure, ” Proceedings of the Symposium | on Earthquake 
and Blast Effects on Structures, Los Angeles, 1952, Published by Earth- 


‘Ss. Levy and Ww. D. ‘Kroll, Introduced by Step- by-Step Integration 
Dynamic Response,” National Bureau of Standards Report, February 


3. Fox and E. Goo oodwin, New Methods for the Inte- 


gration of Ordinary Differential Proc. Camb. Soc., 
Vol. 45, pp. . 373- 388, 1949. 


Timoshenko, “Vibration in Engineering,” D. Van Nostrand 
- Co., New York, first edition, 1928, pp. 79-81. . (and in later editions, see 


4th edition, 1954, pp. 143- 145. 
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ENGINEE RING MECHANICS DIVISION 
ind re- Proceedings of the Societ Civil E ngineer 
only 
ULTIMATE STRENGTH CRITERIA FOR REINFORCED CONCRETE 
greater for ultimate “members are derived by de- 
. These - termining analytically the value of extreme compression edge ‘strain 1 which : re- 2 a 
ne rela-j "sults in a maximum value of moment or load. The following rectangular a a 
on that _ members are considered: (a) homogeneous (b) ‘reinforced concrete _ 
es de- eaten and (c) eccentrically loaded reinforced concrete columns. In all aece = 
of the “cases, ultimate strength so derived is in agreement with tests, and is a se 


function of the stress at the compression and — 


INTRODUCTION 
- Ultimate strength of stuestarel members is ‘commonly determined on the 
_ basis of the assumptions that plane sections remain plane during bending and 
: that stress is a function of strain only. Within the range of linear stress- _ 
strain relationship the stress distribution throughout a member is determined 

_ from the equations of equilibrium of forces and of moments together with the 


¥ of linear distribution of strains. if the stress- ~strain relationshi 


stress distribution, are commonly introduced, or "else the ultimate 
Strength design equations are derived empirically from test data. 


a . By another approach, , the problem | of ultimate strength may be studied ana- 
+ lytically by finding the maximum value of certain load functions expressed in ae { 
terms of the internal resisting forces of the loaded member. . This can be done a i 
"without defining mathematically the stress-strain relationship of the inelastic = 
material. In this manner, criteria for ultimate strength may be derived which © 
_ Discussion open until December 1, 1959. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 2095 is 
« part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings ss 
of the American Society of Civil Engineers, Vol. 85, No. EM 3, July, 1959. 
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are independent of the quantitative definition of the stress-strain 1 relationship 


Considerable difficulty, due to complexities of the mathematics involved, 


of 


may be encountered in maximizing the load function in 1 the case of members — 
_with an irregular cross- section or of non-homogeneous members such as rein- 
forced concrete. However, for homogeneous rectangular b beams, and for rein- 


forced concrete rectangular members governed by tension, this procedure 


leads to relationships | suitable for use in design without simplification. 


_ The notation of the American Concrete Institute Building =r 318-56) 


distance from eitreme compressive fiber to neutral axi 
oS distance from extreme compressive fiber centroid of tensile rein- 


d’ aa distance ee extreme compressive fiber t to. centroid of compressive 


= = 28-day strength of concrete 
= = stress in tensile reinforcement 

stress in compressive 1 th 


= = yield point of tensile 4 


point of compressive reinforcement 


fy = = Few 


distance extreme compete fiber and resultant of 


Itimate resisting moment 


axialload 


— 
— 4h 
~ i 
DB = width Ola section measured raliel to the neutral axis 
= 
||. 


si 
= distance neutral axis 


= strain in compressive ace! 
= strain in tensile fiber 


= in n extreme « compressive fiber at ultimate 


Fundamental Relationships 


Pure » bending of { an inelastic material is considered, for which ‘the Mae 
t is a function of the strain, €, given by f= F(e). It is assumed that stress is 
a function of ares only, and that the stress function is the same for tension a 
and compression, i.e. F(-€) = -F(e). it is further assumed that plane sections 


plastic modulus 


= a member with a symmetrical cross section such as that shown in Fig. 1 
is subjected | to pure bending i in the plane of symmetry, the two equations of - 7 
equilibrium and the equation expressing of are: 


neutral axis, b = Bly). Substituting into (2), 
by Eq. (3), y = and de, we obtain: 


The moment M is thus as a function of one independent 


the material as in flewrel exhibits 
“stress- ~strain relationship of the type shown by Fig. 2 characterized by a a de- 


; “with a corresponding increase in €,, the compressive and tensile forces can- “= 
' not increase indefinitely, but, in the absence of a prior failure, will reach jae 
certain maximum values and thereafter decrease. Also, within the non- -linear os 


“Ta 


— 
| 
nship 
— 
co 
— 
—— 
| 
] | — 
— 
ant of 
itral 
—— 
oh _ of the resisting forces closer to the neutral axis which results in a shortening ie — 
|) 


aditic 
cres 
velc 


shi 
||-STRAIN AND STRESS DISTRIBUTIO ON i: 


f. 


—— 


i 
| 
L STRES S-STRAIN 
| 


CONCRETE 

of the internal moment arm. It is evident, then, that as M increases a con-_ 
dition may be reached beyond which a . further increase of Eo will cause a de- 
crease of the internal resisting moment. At this point the member has de-— wen 
veloped its ultimate flexural strength, My, as is shown in Fig. 3. An analytical _ 
study of these ‘maximum conditions discloses certain relationships which ap- > 
pear to be in agreement with experimental data and which confirm certain de- 
sign a which have been arrived at earlier by more or less empirical — - 


‘The maximum moment i is then a function of the s strain by substitut- 
ing Eq = €y into Eqs. (1a) and (2a). 


‘Rectangular Cross Section 


Ina beam of rectangular cross 


The value of the strain * sala this maximum criterion is designated ‘i 
the strain at ultimate strength. Multiplying E Eq. (5) by and 
= €y and = = fy, | leads to: 


| 
ede de am. 


By ees with Eq. (2c) it is seen that the second term 
= 


The stress fh at the extreme edge depends only on the stress function i ; 
when the ¢ cross section i is It is Eq. (5a), which 


| 
q 
| 
i 
—— 
| le width b= B(y)isaconstan. f 
acts in tne plane of symmetry of the cross section, __ 
‘ 
— 
© 
ia 
‘Sp, of a rectangular cros plastic modulus, 
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Fig PCA ECCENTRIC LOAD 
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ot 
— 
: 


This equation that at maximum moment, , the moment with respect 


to the stress axis of the area under the stress-strain diagram between the | 


given limits is equal to the moment about the stress axis of a rectangle with 


a _ Hence, for a material with an identical stress- strain relationship in tension 
and compression n of the type ‘shown in 1 Fig. 2, the ultimate moment for sym- — 


; metrical bending of a beam with a rectangular cross section may be de- 
termined by assuming a rectangular stress distribution with a stress equal to 
the actual stress at the extreme edge of the cross section. _ This relationship _ 


is true as long as the maximum moment satisfies Eq. (4), asisshownin 


Dividing | Eq. (7) by the stress integral, we obtain: 


a 


left- hand side of tte) is the divided the area of the 
stress-strain diagram, and is thus equal to the distance from the stress axis Ft 7 
to the centroid of the area under the diagram between the given limits. 2 nl 


S 


~ 


, abe this distance, in accordance with common beam notation, as (1 - - kaeu, 


m 


The stress as given | by Eq. (8) for the maximum moment tis a 
ove € k 
which is to fy only v when ke = = 2, , Menon, ‘the rectangular block 
- from which the ultimate moment may be determined is equivalent to the actual e! 
stress block only insofar as it gives an identical value for the bending moment. io 4 
It cannot be used to compute exact values of the total compressive and tensile ¢ ‘ 
forces acting on the cross section, since the stress at the extreme ne edge is 3 not 
necessarily equal to the average stress onthe cross section. 


The total compressive and — — wae on the cross section of the 


“beam are given by the function: 

4 


Eq. (10) € 2/oe and yields for maximum 


| 
arg 
4 
— 
— 
| 
a 


wn 


a ‘Hence, at the maximum value of the internal resisting forces ina member | 

of the type under discussion, the average stress between the neutral axis and — 
the extreme fiber is equal to the stress at the extreme fiber. As indicated by 
- Eq. (9) this condition corresponds to the maximum moment only if k2 = 1/2... “2 

if k9 < 1/2 at the ultimate strength, then 1 fy will be greater than faye, e 


maximum n moment will be reached before C and - 7. their maximum Pat 


The stress-strain ‘relationship of concrete in flexure, used later in the 
Co analysis of reinforced concrete members, was determined at the PCA labo- 
j ratories At a special testing technique developed specifically for that 
A plain concrete specimen of rectangular cross section, Fig. 4, 


purpose. ( 


was loaded by two eccentric loads, the major thrust Pj being applied bya __ 
testing machine. The minor thrust Po provided by a hydraulic jack was varie 
independently in such a manner that the neutral axis was maintained at one a 
face of the specimen. The two loads and the strains in the extreme fibers of 
the specimen were measured continuously throughout the loading range to 4 


failure, so that the average compressive force on the cross-sectional area, 

as well as the position of the resultant of the compressive forces, could be — 
_ correlated with the corresponding strain in the extreme compressive fibers. 
_ The stress-strain relationship was then derived by mathematical analysis. a 
- = The conditions in the specimen described above can be compared to those 


= relationship in both compression and tension identical to the areme- 9 
of the concrete in f flexural compression, asis shown in 

‘It then follows that the moment of the compressive forces in the specimen 
about the neutral axis is equal to one half the moment in the hypothetical rae 
_- if the curvature of both is the same. | Accordingly, | the maximum 


get: 

_ From the test data of the twenty specimens tested at the PCA laboratories, 
_ and the stress-strain relationships derived from them, the values of the maxi- 
mum moments and corresponding stresses in the extreme compressive fibers 

i were obtained. Substituting these stress values into Eq. (12), the theoretical | 


values of the maximum moments were calculated. The results, 


stress- strain in relationships in ( 


Pe In a rectangular reinforced concrete beam, Fig. 6 - te Sal the stress in | 


au fensile and reinforcement is equal to the stress, the | 


3 
| 
a > 
a 
| 
_§ 
“ 


fy 
Fig, CONDITIONS AT ULTIMATE LOAD 


‘PCA ECCENTRIC LOAD 
> 


Axis” 


| 
6- - CONDITIONS IN R.C. 
YIELDING OF REINFORCEMENT 


_ BD CONCRETE 
— 
ne a fi 
nose q — 
— 
fibers 
— 
— 


+ ‘TABLE 


ULTIMATE MOMENT IN P PCA ECCENTRIC LOAD ‘SPECIMEN 
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‘The moment arm of the resultant of the ‘concrete compressive forces with | 


respect to the centroid of the tensile r reinforcement isdeterminedto be: 


F(e)de - F(ejede | 

—_ ing Eq. (15) into Eq. (14), differentiating with respect to €,, and equat- 


‘ot 


Dividing Eq. (16) by “ Fle)de | vields 


and 
p- -p = + 
My 


__‘The stress f, in the extreme concrete fiber at ultimate strength is aS, &§ 


a of the stress-strain relationship of concrete in flexure only, when the 


beam cross section is rectangular and when the beam is governed by tension. ie fl 


‘The stress fu be from (17), which may be rewritten inthe 


i 
q 
est 
Ay 
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J uly, 


= 2k 


’ - ‘Hence, when the maximum moment is reached in a rectangular reinforced | q 
Pod concrete beam governed by tension under symmetrical bending, the concrete — 

; stress distribution in the compression zone is such that the total compressive 
force is - equal to that obtained from an equivalent rectangular stress block in ¢ 
which: (a) the stress is equal to the actual stress in the extreme compressive 
fibers, and (b) the depth is equal to twice the distance from the extreme com- © 


pressive fibers to the centroid of the compressive forces. . Thus, the 


ne: For ‘design purposes, it is then sufficient to establish ‘the lt of the stress — 
fy as a function of concrete strength fo. This can be done by solving Eq. (18) - 


. ‘ when test values of beam properties, concrete strength, and ultimate moment — 


‘ae known, in which case it is not necessary to know the entire concrete 
From Eq. (19), the stress and the corresponding strain can 
termined when the stress- -strain relationship of concrete is from ex- 
‘periments. Using the stress-strain relationship obtained by Hognestad, 
4 Hanson and ‘McHenry, (1) a relationship between fy and f¢ was determined for 
concrete at age 28 days. The function obtained is represented by the curve 


shown in Fig. 7, ont may be very by the equation 


if. 
in which fa and are in psi. 


7 _ which is identical to, and therefore confirms, Eq. (A3) of the 1956 ACI Building 


Letting Ag = 0, and substituting Ag = pbd and q = Ply/t Ba. (22) yields 


Tension Failure of Becentrically Loaded Short Columns 


4 
— 
beam moment given b Eq. (18) is not sensitive 
to the value of | ndesign. A straight- 
strength and is approximately 0.0023. 
is denned herein as a Column Of sien proportions, and with 
Yk... wet ‘a an axial load applied at such eccentricity, that the eccentricity may be con- ff 
If such a column of rectangular cross section, shown in Fig. 8, is loaded 
— 


ake 


Concrete 


Vee 


Stress in extreme edge, f,, in ksi 


7- “STRESS IN EXTREME EDGE AT ULTIMATE 


STRENGTH OF R. C. MEMBERS" 


| Assuming that both tensile and compressive reinforcement yield before 
sive failure, f, = fy and = fy. The two equations of oad 


y 


Figs 8- CONDITIONS 


LOADED R. Cc. COLUMN 


wa 


— | 
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“The moment arm of the of the co concrete compressive 


= 
Solving Eq. (23a) for c and id substituting 1 into nto EG. (28) w 


in: 
(20) with respect to equating to zero, and simelifying 
/ Hence the criterion derived te: maximum moment in rectangular rein- — 
forced concrete beams failing in tension applies also as a criterion for the 
maximum eccentric load on rectangular reinforced concrete short columns 
ta By substituting Eq. (17) into Eq. (29) the length of the moment arm (d - ke) | 


at ultimate strength is determined. With this value of (d - ko") the ultimate va 
strength P,, is then obtained by solving Eq. (24a). As p'bd, 
and at {, by means of Eq. (21), this leads to: 


which is equivalent t to, ‘and ther contirms, ‘Eq. ‘as) of the 1956 ACI 


Building Code. i 

In the case of materials with a pelationshie exhibiting 
scending branch following the maximum stress, the ultimate strength of _ 


= members can be determined analytically by maximizing the en 


able strength function. In this criteria for ultimate strength are 


o The maximum resisting 1 moment of rectangular eanapenceus' beams in 
rroduct of the — of the 


— 

i 

— 
— 
— 

a — 
q 


EM 
and the stress in the fibers at ‘ultimate If the stress- 
: strain relationship of the material is known, the stress | fy in the extreme [| 
* fibers at ultimate strength may be found, since the moment with respect — 
to the stress axis of the area under the stress-strain curve at ultimate > 
strength must be equal to the moment of a rectangle with the ultimate 3 
strain €y and fy = F(ey) as its 
(2) The ultimate strength of rectangular reinforced concrete members | 
_ governed by tension isa function o of the properties of the cross section 
of the member, the yield strength of the reinforcement, and the concrete 
stress distribution as represented solely by the stress in the extreme _ 
_ compressive fibers at ultimate strength. This function is expressed a 
(22) for beams in pure bending, and by Eq. (30) for. 
loaded short columns. The concrete stress in the extreme fibers at 
ultimate strength is euch that the area under the corresponding ool 
strain curve is equal to the area of a rectangle with fy = Fleu) and 
(3) The relationship between stress in the extreme compressive fibers at — 
on ultimate strength, fy, and the 28-day concrete cylinder strength was ; 
found to be approximated with | sufficient accuracy by f, = 0.85f¢, corre- 
sponding to a strain of 0.0025. The analytical criteria presented herein 
then lead to design equations identical to those given by the 1956 ACI 
_ Building Code for ultimate strength, governed by tension, of members © 
_ of rectangular cross section. Hence, the study presented reflects a 


- theoretical substantiation of Eqs. (a3) and (A9) of the ACI Building a 


- The paper presented herein was submitted by the author to the | Graduate — 
- School of Northwestern University as a project report in partial fulfill ment of 
the requirements for a Master of Science Degree in Civil Engineering. = 
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» SWELLING AND CREEP IN CEMENT 


_ A. Hrennikoff,! F. ASCE 


paper has been prompted by an earlier study of creep. (1) It has been observed 
that both phenomena are associated with the presence of water, and that they yy 
both manifest themselves in linear and volumetric changes. This suggests ‘ 
_ common origin for the two phenomena related to water. "The shrinkage study 
_ has led to the enunciation of the active water theory, which is described in 
ar paper in detail. The theory attributes the phenomena considered to the 
stress condition in the film of water surrounding cement grains when cement 
is moist. | ‘The theory agrees well with several aspects of cement behaviour in 


It has long been known that shrinks on swells on wetting. 
_A detailed study of the circumstances surrounding : a cycling variation of ing. 
content in cement provides some clues on the role of moisture in cement per- | 
mitting to link the phenomenon of shrinkage and swelling to creep. i 
a The experimental study of shrinkage and swelling involved preparation on of ; 
_ prismatic specimens of cement approximately 3/4" x 3/4" x 4" in size, of 
different water cement ratios and cement grain | sizes. . These specimens were 
. subjected to cycles of ‘drying and wetting with water and other liquids. ‘The 
| principal items of information obtained from the iemneuae were the values” oy 
of weight and length | taken periodically in the test. 
_ The weights were taken at first by an ordinary ay to 0.01 g gm., and es 
later by a self indicating Mettler balance to the same precision. _ The length 
_ Note: Discussion open until December 1, 1959. To extend the closing date one 3 
a written request must be filed with the Executive Secretary, . ASCE. Paper 2096 is 
_ part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings ~ sas 
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was ‘measured directly by a home made apparatus, shown in the photograph, ye. 


“Fig. 2 involving a frame and a dial gauge, sensitive to 0.0001". Although - 
somewhat crude, the method of the length measurement was found Yn 
A surprising degree of consistency was observed in successive measurements 
of length, although the individual readings could not be relied on closer than 
to 0. 0002" or 0.0003". The errors were caused mostly by the relative rough- 
ness of the ends of the specimens and by the difficulty of bringing the measur- 
ing points of the apparatus exactly on the targets 
For interpretation of the results, the readings of the increments or dec re- 
ments of the weight and of the length were plotted on the graph paper asthe _ 
abscissae and the ordinates respectively. The increments were plotted in the © 
‘positive directions of the axes, i.e. to the right and up, and the decrements in | 


Each graph refers to one specimen only, and it represents its several 
cycles of wetting and drying. The loss and gain of weight is looked upon as 
the decrement and increment of the internal volume of voids occupied by the 

liquid. When the liquid is water its specific gravity is assumed unity, although 
it may not always be so, and the losses and gains in weight are plotted as the - 
Variations of the internal volume of water incu.cms. 
Other liquids: kerosene, lubricating oil SAE.10 and methyl alcohol were 
also used. The changes in weight incurred by these liquids were converted _ 

into volumes, ‘prior to plotting them on by their specific 


Several types of cement "studied: ordinary cement, fine 


the negative directions. 


cement, coarse cement, of water-cement ratios corresponding to plastic and 


4 
more liquid consistencies. Several specimens of each type were used, and 
were subjected to different manners of drying and wetting. Some specimens | 

were hardened in the moist closet under spray, others in the atmosphere of 


i 


: steam in the boiler, under conditions similar to the ones prescribed by. ASTM 
‘Figs. 2 to ‘T represent a set of tumheeh graphs of this kind, the proportion o of 
_ ingredients by weight being 0.003 : 0.24 : 1 of pozzolith (a dispersing agent), a 

water and ordinary portland Type 1, manufactured by British 


, closet, removed from the forms and then kept for two days (about 6 hours Ra 
each day) in the atmosphere of steam. This was followed by drying them in 
4 “y the oven at 120° F, and 130° F. When ina week or 10 days the weight was sub- 
stantially stabilized the specimens were transferred to the oven with 
a temperature of 200° F., and kept there until the weight mane to decrease. — 
‘This followed by wetting in several different nt ways. 


we The weights and lengths of the prisms were taken at appropriate time inter- | 


vals, more frequently when they were changing fast, and the ages of prisms 
a in n days at the time of measurement were recorded on the graphs. When dry-— 


_ Here are some of the significant features of the graphs: of 


increase in weight. The latter is caused by hydration roe in decrease 
of the absolute volume of products of cement hydration, followed by a tendency 


to form voids and by the suction of the moisture from the outside. This phe- a | 


= - nomenon, not particularly important in connection with the present studies _ 
by others and described by the term “self dessication” (2) 


os ing in the oven the e length m measurements were naturally taken after ‘cooling the : 
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Drying. In spite of minor the curves 
the curves is their flatness at the early stages of drying and steepness at = 
later stages, indicating that the loss of water at the end of drying results in a 

greater shrinkage per gram or C.c. of water lost, than at the early stage. Ex- f 
; _ pressed differently, the water lost early in drying is comparatively “inactive” ay . 
as far as il is concerned, whiie the water lost later is “active” a 


4 all six specimens are very similar. The most significant characteristic of 


Structure of Hardened C Cement 
mostly of gel, (3) a hard sponge like mass, 
whose cells or individual units are too small to be seen under an ordinary — 
if 


microscope. Gel contains some 20- 25% of voids filled mostly, with water, 

cement has been kept moist. 
Hardened cement also contains some ‘unhydrated material—remains of 
larger cement grains, and some bubbles of air trapped in cement during its 4 
making. The air bubbles contain no gel, because gel is a of alienate a ar 


a. Since gel is amorphous it should be visualized as a of the elements a 
of hydrated cement surrounded by water according to the scheme of Fig. 8 


| and not in the form of a solid mass containing intercommunicating internal = 


_ cavities according to Fig. 9. This follows from the fact that gel is deposited = : 
in water and in view of being amorphous, its elements do not merge into one ‘ 


another but remain distinct, and consequently separated from each other by a 
layer of water, however thin it maybe. 
‘The arrangement of Fig. 9 could be true when applied to a single crystal. Lo 
> The voids then would be the imperfections in the crystal brought about by some — 
disturbing influences. . This however does not apply to cement. Some . 
vations with electron microscope(4) corroborate the structure of Fig. 8. 
sl Considering the graphs of drying in conjunction with Fig. 8, ‘it appears | at 
- the inactive water belongs to the enlarged parts of the gel voids, because 
_ these parts would dry first, while the water coming from the ——- of grain 7 
to grain contacts belongs to the active kind. 


“After aie, te specimen 7 § 3 (Fig. 4) was wetted by immersing in 
water. By the next day the specimen had absorbed nearly the full weight of ss 
- water lost earlier in drying, but had recovered only a small part of its shrunk — 

— length. The curve of the first day of wetting is accordingly flat. Continued — 
immersion resulted in a slow gain in weight, accompanied by a substantial - 

‘Thus the water absorbed early in the course of wetting was ~ 
while the water absorbed later—active. it appears then that the swelling is ay 
_largely due to the saturation of the less accessible voids in the gel, , i.e. the 7? 

of the grain to contacts. 


specimen 78 5 (Fig. 6) 
sulted in absorption of moisture. The wetting of the specimen has been stow 
and accordingly, the | curve of wetting in air is quite different from the curve 
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somewhat to the right of it. means that, apart from a amount 
inactive water, the water absorbed first was active and the water absorbed _ 
later—less so. The moisture absorbed early thus found its way into the 
_certain amount of moisture accompanied by recovery of half of the ——— lost 
drying, further absorption of moisture from the air nearly stopped. 
‘The specimen 7 S 4 (Fig. 5) was moistened in a closed can containing some 
a water at the bottom, below the specimen. The can was kept in an oven at a 

140° saturation in this case was naturally fast because the hot water 


vapour was than the cool ve vapour. The = the curve 


with Inert Liquids 


_ they normally do not affect the Aint and other mechanical properties of “7 
_ cement to any appreciable extent. (5) Accordingly their lines of saturation on 
= the graphs are almost completely horizontal, i.e. these liquids cause no swell- 
ing of cement, although the internal volumes comparable 
to the volumes vacated by water on drying. E. ee 
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After saturation with these liquids, which required a considerable 
of time especially in the case of oil, the ; specimens were placed in water and E 

7 here a surprising phenomenon was observed: the specimens began to some p 


water: without displacement of kerosene or oil, and at the same time they be- 
gan to swell along a comparatively straight line more or less parallel to the _ 
line of drying at 200° F. In due time the specimens swelled beyond their origi- — 
’ nal! length before drying, and absorbed a greater total total volume of liquid thar than wa ; i 
‘This experiment with successive saturation of hardened cement first with 
== liquid and then with water is the most significant in the whole series of 3 
tests and it gives a clue to the nature of active water. How can water pene- 7 
. _ trate cement, if all pores in it are filled with kerosene, a liquid nonmiscible r 
_ with water, and especially when no kerosene is forced out? An explanation fi. ’ 
. suggests itself, that the water enters the cement not by the way of its voids ir. ‘ ‘ 
- but by the way of creeping along the surface of the cement elements. On the =; 
4 “surface of ‘the specimens there are some ‘Spots, from which the ‘kerosene has 


grains. Once placed in contact with some cement grains, the water creeps _ 

around them coming in contact with more grains where the latter touch ir? 
other, and in this manner penetrating inside the cement mass. ete 
___- Virtual equality of the moduli of elasticity and of the Poisson’s ratios in 
Bs the moist and dry cement specimens, demonstrated in the previous studies, (1) 
5 suggests that the presence of water does not soften the grains of cement gel 7 
or change the grains themselves in any other way. On first thought these ete 

a _ seem to conflict with the phenomenon of swelling, but the hypothesis ofthe 
a ; water creeping along the grains and spreading them apart like a wedge, while 
; _ jumping from one grain to another, provides a plausible explanation. Accord: 


molecular forces transform this water from liquid into solid—a plastic or i a 
he 4 viscous solid rather than a rigid solid, since a rigid solid could not creep 
along the grains. The water that fills the central parts of the voids in Fig. 8 

is hot active because it is too far removed from cement grains | to be: attracted 


. ie inert liquid makes no difference with regard to swelling of the cemen 

as Judging by the increase in 1 weight in of swelling in Fig. 2, 
weight of active water equals approximately one third of the weight of the total 
_ chemically uncombined water. . Other tests make this fraction even smaller. a 
_ Separation of cement grains by the film water occurs against a strong inter- 
| attraction holding them together. The active water then must be 
under a very strong compression. . For this reason active water is entirely Te 


| = from capillary water which is in the state of tension or, more iw i 


correctly, in the state of compression smaller than atmospheric pressure. (6) a i 
Consequently, regions occupied by active water similarly to the ones occupied — 
y the cement grains themselves, are unavailable for « capillary flow. Only the 
parts of the void spaces filled by inactive water determine the permeability — 2 

_ of cement mass to water flow, while full, although diminished by shrinkage, 
void spaces are available for kerosene flow. 
conclusion that the volumetric changes in cement are caused by water 

Present ina state other then conflicts with the 


4 
4 

| 
to UMS View, active water 1s the Water lorming 4a tim around the ceme 
grains or adsorbed by their surface. This water is not an ordinary liquid ij a 
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_ _ The writer experimented with placement of small drops of water and kero- 
_ -Sene on | the surface ce of hardened cement, ' The surface had been prepared by ie a 
: drying a specimen, cutting it and polishing roughly the surface of the cut. This 
surface was held in a horizontal position and viewed through a microscope = 3 

pointed The kerosene was: observed to ‘spread quickly over | Wh 


area or to contract to a ‘smaller area. At the same time the water in the Pet 1) - 

percolated into the cement and spread slowly underneath the surface. As the eal 

_ drop decreased in volume the inclination of its sides remained the same but | . Ki mt 

its skin formed wrinkles, being apparently supported by a crystalline frame-_ 

precipitated from the solution. When the drop finally it left ak 

‘The point which has a considerable importance in appraising the roleof | _—_the 
Pi. capillary forces in moist cement is the magnitude of the angle of contact of ta : 

the water surface with the cement surface, which was found to be 80° - 90°. jis: 7 

This compares with the angle of 15° - 20° formed by the water in similar a 
circumstances with the ee surface. In the textbooks the latter angle is 

usually assumed zero. ( 7) The contact angle of the liquid has a direct relation _ 

_ to the height of its capillary rise in small bore tubes and to the negative water ay 

_ pressure formed underneath the surface meniscus. If in a fine glass tube 3 $ ‘int 


capillary rise is high, in a tube of the same diameter but made of cement it 
is smaller by a coefficient equal to the cosine of the contact angle. (8) 8) Should 
_ the water in the fine passages between the elements of cement gel meet the 
cement at the same 80° - 90° angle as at the outside cement surface, the 
- capillary phenomena associated with the ; presence of water inside ‘cement 
would be insignificant. On the other hand the kerosene, forming a zero angle 
= with the cement surface, should produce a substantial capillary rise in cement 
q - , even though the surface tension of kerosene is only half as great as that of a 
3 Experiments were made with the rise of kerosene and water in dry cement 
prisms whose bottom ends were immersed 1/16" into these liquids. The rise — 
of kerosene was the faster of the two and it could have been caused by the -_ 
capillary forces only; it produced no contraction or expansion of the prism. 
The water rise, . however, although slower and smaller than that of the kero = 
_ gene, was accompanied by a substantial swelling. The kerosene rise proves ¥ 
_ that swelling of cement is not caused by capillary forces, while swelling con- 
"current with the water rise suggests that such rise is brought about at least 
. in part by a cause other than capillary. This cause, accounting a also forthe 
The presence of film water in was demonstrated by several other 
t investigators, (3) without however attributing to it the mechanical action = 


on ct 


= 
Be 


causes of the phenomena discussed. 


22 sthe cause of swellingand | 
illary action of water is 
holds that the capillary 
"| 
| 
: 
| 
| 
| 
a 
] 
|, 
| 
| 
| 


nt 


q Water 


a a sufficiently large positive number. The law of attraction of this type | 
the force decrease rapidly as the distance increases. 
_ In addition to the attraction there must also be a force of repulsion R, be 
cause once the particles come in contact they are quickly stopped in their 
- mutual approach. Each particle is then in equilibrium under the forces A and | 
i The > expression for the force R may be assumed similar to Eq. (1) except nan 
that the exponent n in the denominator is even greater for R than for A, since i 
the force of repulsion appears to be virtually non- existant a 
Presenting. the force- distance relations graphically in Fig. the curve 


“The distance dp at which the particles find their is by 


“ternal forces F, shown dotted in Fig. 10. The total force tending to bring the 
. ‘particles together is now (A+F), while the force resisting an unduly close ap- be 
proach of the particles is still R. . if the ordinates of the line A are increased . 
_ by the force F, the new distance of equilibrium dy (Fig. 11) is found by the = 


“ultimate tensile strength for one pair of particles. 3 


| This general qualitative scheme of mechanical action of dry fini 
a under stress can scarcely be doubted. Now the water will be intro- 


In Fig. 12(a), AB represents the surface of a cement grain in contact with — 
water, pictured plane for simplicity. It is assumed that the water is — 
| racked by cement in accordance with an inverse power law of the type of | 
Eq. (1). ‘The variation of this force with the distance t from the cement — a 
surface is pictured in Fig. 12(b). This force f, expressed in units of force ee 


unit volume of liquid is very great near the surface AB and is virtually zero | 


_ Attraction of water to cement creates pressure p in the former. Itsin- _ fe } 
i 

can be easily determined from equilibrium of a prism of water i 

CDLE in Fig. 12(a). prism of cross-section dA extends fromt tothe _ } 


at some distance to. It is only within the thickness / that the liquid is active. <a | 


er aa grains like other elements or particles composing solid bodies are ace 
held together by intergranular attraction. The law of attraction is not known 
| 
i 
same time the graph of R must become nearly horizontal at its bottom end, 
4 
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he 
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In order to be mathematically exact the e upper limit of the integral s should a 


The variation of p in accordance with Eq. (2) is shown in Fig. 12(e). 
curve of p resembles the curve of f in its general shape. 

It is evident from the expression (2) that the pressure OEE represent 

_ the areas under the curve of the force of attraction. Thus the pressure p at 

the distance t equals the area KTU and the pressure Pi | - the: area K Ti — 
a : a It has been assumed so far that there is an unlimited amount of water avail- 

ae in contact with cement, in other words the thickness of the water layer = 
is ; not less than to- Suppose now that this thickness of water is t < i... This — 
_ does not change the graph of f, although its effective part is now only U US, x 
but the pressure graph change _ 

-. at the distance t ty it becom sored 

In dther words the | pressure in Fig. -12(c) a are n now ww mannered from. 
ia What has been said about the attraction of water to cement grain with a i 
"plane surface is believed to apply without a change to a curved cement surface — 
eben that the radius of curvature of the grains is much greater than the a 
_ thickness of layers of adsorbed water. If however the curvature of the grain — 
ts aly eye to the thickness of the film, the situation is probably different. — 
- Once water comes in contact with cement and forms a film, it creeps along, 


“ness. ;. The cause of lateral ‘spreading s of the film around the | grain is the need 7 


_ slowness of swelling of dry cement prisms on immersion in water, peeementne-— 
; _ ly after saturation with inert liquid. In this respect film water behaves as —_ 7 a 


a for lateral confinement. Thus, the water prism C DL E considered above in = 
si in Fig. 12(a), being under a strong pressure in horizontal direction, is not in — 
a : equilibrium unless it is restrained laterally i. -€. in up and down direction, a 
’. similar neighbouring prisms under pressure. In the absence of such con- ae 

a liquid of having its pressure at a point equal in all directions and normal to 
‘7 the plane of action. Film water differs however from ordinary water in that 
” :. respect that a period of time, and even a long period, is required for _ 
adi _ justment of its shape to a change in conditions. This is indicated by the biey ai 


plastic solid. This change in the properties of water should be attributed to 
strong forces of attraction bringing about some regularity in the arrangement is 
Let the curve in Fig. 13 represent the diagram of variation of water 
7 pressure in the film around thé cement grain, similar to Fig. ha a 
Ed pressure Po next to the surface of the grain may be relatively high. Divide Po . 
_ into a number of equal parts Ap. Four such parts are shown in Fig. 13. The a 
 abscissae of the points with the ordinates 4Ap, 3Ap, etc., are tp = 0, ty, | 
_ Imagine now an isolated spherical grain of cement G surrounded by ample “i 
mount of w water (Fig. 14). The pressure in the film water of the g grain may — 


a 
| 
— 
— 
— 
| 


appear in as circular  isobars. “The is isobar zero with 

_ pressure 4Ap, is the surface of the grain, the next isobar ni number one has the 
pressure 3Ap, and the distance t, from the grain. The distances between the 
_two successive increase outward. distance tg _4 between the iso- 


- grain is limited to the thickness to for example, the isobars number 3 and i 
disappear, while all others remain in exactly the same positions, although — 
their pressures decrease by 2Ap, i.e. the surface of the grain becomes call 

eS 2Ap, the isobar #1 becomes Ap, and the isobar #2 becomes zero. 

sdf a plane, such as AB (Fig . 14) is passed through the film water the Fit 
a _ pressure in this plane is normal to it, and the variation of its intensity is ai 


ah. Suppose that an inert | body Tis brought into the vicinity of the grain a 

(Fig. 14). The isobars outside the region I including the space between I rm 
G remain completely unchanged. The reason for this condition is that although 

the thickness of the water film between G and I is reduced, the water on the — e 

; j ; sides ¢ of I crowds in and maintains the same pressure between I and G as else 
where at the corresponding distances from the grain surface. 

q : 7 ob Le If two cement grains § surrounded by skins of film water are brought in cheat | 

proximity within the range of their intergranular attraction, they approach — 
each other and find a position of equilibrium by distorting the films at the 
point of their near contact. ‘The’ system of isobars in that vicinity is complete- 


a ly modified, and the water pressure in the film transmits not only the ee. 7 
traction between the cement ; and wi water b but also the attraction between cement 
interaction of dry | grains of cement discussed earlier the force of re ke 
pulsion of the grains R (Fig. 11) is one of the two principal factors. As to the 
moist grains, considered now, the repulsion force is not believed to be of ‘<a 
primary significance. If the grains are smooth the water film islikelyto 
prevent their touching each other, which is the prerequisite for bringing the 
force R into action. Should the ) grain surface be jagged it is possible that __ 
- some solid to solid contact may take place at the protuberances, while the “a 
major parts of the film remain unpierced, however the action of R under these a 
Iwo se is probably of secondary importance, and may now be ignored. umes 
_ Two sets of isobars may be considered when two cement grains are brought 
M 4 in close proximity: the isobars corresponding to the independent action of ie 


— the produced by the combination of attractions of both 


4 se R i concentric circles, and so they remain when one grain is brought close | ie 
onl ‘another, as has been shown in Fig. 14 by reference to the body I. 
ss The two sets of isobars may also be combined into one by direct addition — 


_ of pressures brought about by the attractions of the two grains. . That this pro- 

_ cedure is correct becomes clear from the following two considerations: first, 
that the forces of attraction of any element of water by the two grains are en- — 
tirely independent of each other, and second, that the pressures produced at nf 

any particular point by each set of forces are equal on all planes, as well as 5 

eS ee _ How this combination of isobars can be carried out is shown in . Fig. 15. 
9 Here three circular isobars are shown surrounding each cement grain. They © 
to pressures Ap, 2 2 (4p) ‘The isobars 4(Ap) 
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CEMENT 
ith the po es of th the grains. Parts of the isobars 0.1(Ap) are also shown. 


The ‘combined pressures can be easily determined at the points of intersection 
_ of the isobars of the two families. Thus, at the intersection ofthe isobars _ - 
2(Ap) and Ap the pressure is 3(Ap) etc. On the straight line between the | a 
centres of the grains the pressure varies between approximately 6(Ap) next es. <9 
to the grains and perhaps 5.5(Ap) half way between the grains. The combined ~~ 
-— easily sketched from these points, are shown by solid lines. ll 
Although the combined isobars | present a clear picture of the pressure | : , 
‘distribution in the film water, the individual circular isobars of Fig. 14 are 
more convenient f for the purpose of calculation in | problems involving swelling - ac 


and shrinkage « of cement, and for this reason they will be use 
e- The oe puatiens in shrinkage is finding how far apart the cement grains 
-would move when seeking their position of equilibrium in conditions of a varia- a 
_ ble thickness of water film or a variable moisture content. . When the water 
film is thin the grains are relatively : close » together. _ As the thickness of the oa 
film increases the grains move away from each other. . The amount of ac 7 


* The exact numerical solution of this problem at present is <a due 
_ to the lack of knowledge of the exact laws of attraction between cement and — 


water, and ‘cement and cement, but a detailed outline of the ‘method follows 


Clearly from the theory presented. 


af Suppose that two equal cement grains (Fig. 16) surrounded by water are 
) into their ‘relative position of equilibrium lying within the field of = 


ori G. ¥ They are the attraction hives A developed by the other grain G2 


and the pressure a the surrounding water film. ‘The latter is produced _ 


The attention of water by cement of this grain. 
a b) The attraction of water by cement of the other ae 


was explained above with reference to Fig. 15, or considered separately, 
which is more expedient inthiscase. 


The water pressure on Gy corresponding to the attraction (a) is self 

_ balanced (Fig. 16b), ice. it produces a zero resultant on the gr grain Gy. Then 

the water pressure on the grain Gj caused by the grain Gg is the only one 

which remains to _be*considered, and it is balanced by the intergranular at- 
traction AL This pressure corresponds to the isobars of the grain Gog at the 
7 surface of the grain Gy (Fig. 16a) and it is depicted in Fig. 16(c) as the line ig i 

_ cbabe. The intensity 0 of this pressure is the greatest at the point a, and the ey 


grain is so far from bag oe G2 pressure aR is always 

ay The water preasert acts on the | grain in the direction normal to its surface, 

_as indicated by arrows in Fig. 16(c), but, as in similar conditions inelementa- 
hydrostatics, , the horizontal components of pressure are mutually balanced, 
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The two pressures caused by these two attractions can be either combined as 
rst, 
4 


and the vertical up a be found by sum-_ 
- ming up or integrating the pressure intensities over the horizontal diametral — 
area of the grain. This resultant is represented by the volume of the body of — 
revolution under the surface described by the pressure line in Fig. 16(c). As 
pointed out ut above, this this pressure produced is equal to the the force of 
If the moisture content of cement is dec reased, the films around the e grains ; 
bensene thinner and the ordinates of the pressure line on the grain ests om 
16c) becomes. smaller, ‘so that the resultant repulsion becomes repre- 
sented by a smaller volume, such as the one under the pressure line bab. — ie . 
Since this force is less than the attraction J A, the grains must come closer to- to- 
_ gether. This raises both the attraction A and the pressure : repulsion, the te | 
latter faster than the former. In this manner the equilibrium is re-established 
with the grains brought closer 
opposite effect occurs with the introduction of more water in a. com-_ 
geile dry cement. This thickens the active films, causing an increase in 
_ their pressure, which brings about greater repulsive forces and some dis- | 
"placement of the grains away from each other. The effective intergranular 4 
: attractions extend apparently for greater distances than the pressure produced 
repulsions and so the active water cannot cause a 
or the slaking of cement 


_ As was pointed out earlier, the volume - shrinkage graphs s (Fig. 2- 
are very steep on the dry end, indicating that at the conclusion 0 of ee = war 
; ‘small loss of moisture corresponds to a relatively large linear r shrinkage. _ i 
However even at this stage, if the loss of moisture is translated into a de- 
crease of thickness of the active film, this decrease = exceeds the ie 
by” which the grains / of cement ap approach each other. Stated ‘differently, if in a 
the process of drying the active films around the grains should become thinner 
say by 10-6 inches, the grains would approach each other in the « course 
— not by me inches but several times less. i 


films of thickness t, small in relation to R. At the points of contact ofthe 


of the smallness of the thickness t the part of the surface area of the grain 
covered by the film of reduced thickness may be considered negligible com- 


Take the volume of cement prism as 37 cc. and the volume of uncombined — . 
water in it removed by drying at 200° F., as 9 cc., which figures are typical f[ 


EE the films are considerably thinner than | elsewhere. However, in view 3 


a of all specimens in the shrinkage graphs, Fig. 2 to Fig. 7. Then the omega | 


volume of dry cement grains may be taken as 37-9 = 28 cc. Actually 
"some uncombined water undoubtedly still remains in cement when it is dried : 


prior to this drying the specimen still ‘contains the volume of water 


all of which is active water, the oe = the film t of this water may ” s j 
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CEMENT 
which, ignoring terms with higher oft 


Strictly R in this expression cement grains 
E the thickness of the water film unexpelled by drying at 200° F., but 

in view of the negligible thi — of t the latter, R R may be taken as the radius 


Apply this formula to the specimen with the graph of Fig. 6 whose drying 


this: ‘specimen 1 at 200° F may be assumed complete ; at 18 days. The volumes ; Be 
of moisture ‘still present in it at 15, 13 and 7 days respectively are 0.25, 0.91 
and 2 2. 55 cc. 


_ Then the film thicknesses by Eq. (4) at different a 


sare: 


0.011 R 


al The decrease in the film thicknesses in diffe in different time intervals are a 


tba 


— 


13 to 15 0,008 

(15 to 18 - 0 003 R 


The corresponding magnitudes of shrinkage in the same pon ane are 
- found from the graph as 24(10-4), 32(10-4) and 20(10-4) inches on 4" length 


or 0.0006, 0.0008 and 0.0005 inches per ‘inch respectively. eee on: 
}  -—Had the specimen shrunk for the full amounts of the decreases in the 


thickness of the film the corresponding shrinkage strains would have been 
-019, 0.008 and 0.003 inches per inch. Thus, the of the im de- 


particular sp specimen the film decrease— ratios for different 
intervals : are as follows: 


13 days - 


0.0006 ~ = 31.7, here the is 


ar 45 to 18 days 0005 0.0005 


a There is considerable variation between different specimens with regard mee 
tot the of lines, from errors 
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observation, partly from causes to be p pointed ou some e speci- 
mens the terminal film decrease - shrinkage ratio is less than 6. vee z 
If it were possible to demonstrate on the basis of the active water theory 
that the film decrease - shrinkage ratio should necessarily be of the order ob- 
served in the tests, or at least much greater than one, that would be an | 

_ ditional evidence in favour of the theory. As the matters stand now, the fully 

"conclusive proof of this proposition is at present out of the question, for lack — _ 


| of the necessary quantitative parameters, but the plausibility of it is easily 


demonstrable by choosing suitable data. 
Fig. 18a shows two cement grains of radii R surrounded by films of water — J 
_ of equal thickness on both grains, although only the film belonging to the grain be, : 
‘ G2 is shown. This film is described by the spherical surfaces 1, 2, 3, 4 and 
4 5 which have equal pressure intervals between them. Thus, ifthe water film 
covering the grains extends to the surface 4, the pressures on the surfaces 
4,3, 2and1 are 0, Ap, 2 (Ap) and 3(Ap) respectively. The film thicknesses 
measured from cement to the surfaces 1, 2, 3, 4 and 5 are assumed arbitrarily 
ty = 0.01 R, tg = 0.025 R, tg = 0.55 R, tg = 0.125 R and ts = 0.30 R respective- 
These values are probably too great in relation to the radius, but smaller 
caamadion would be too small to handle in a graphical demonstration. = = 
® _ Suppose that the water film on both grains e extends to the surface 4 and that 
the position of equilibrium of the grains covered by this amount of water 
such that the grain hy touches the surface 1 of the grain G2. . This is the po- 


‘Fig. 18(b). It is constructed by observing the fees) of intersection of the - 
surface of the grain G, with the isobars 4, 3, 2 and 1 of the grain Gog and by <a jie 


measuring pe ae of these points of intersection from the vertical axis 5 a 


_ -The force of repulsion between the grains is equal to the volume of the body i 


Bas revolution, , whose axial cross-section is represented by the pressure dia- 


By the theorem of Pappus, assuming the shape of the generating figure as 


a triangle and scaling its base, the force of repulsion “ found thus: ae 


p 345 R) 2m (1/3)(0.345) R = 0.3738 (4p 

‘Since the relative position of the grains considered here is assumed to be the 

position of equilibrium, the force of attraction A of the grains in this anal 


Assume now that after ofa a amount of water ome grains 

a move away from each other, so that the surface of the grain Gj touches the _ " 

: sphere 2 of the grain Go. _ The question is, how thick a water film is required 
for the equilibrium. Assuming tentatively that the necessary film extends 
the way to sphere 5, the water repulsion then corresponds to the pressure dia 

: gram shown in Fig. 18(c) and extending vertically from the isobar 5 to the ” 
isobar 2. Iti is evident however by inspection that this repulsion is too great _ 

it exceeds the one corresponding to Fig. 18(b), while actually it should hind 
be smaller than it, because the force of attraction between the grains = Be 
greater distance apart must be | smaller. Therefore the increase in the 
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1989 
Revising the earlier assumption, assume the outer boundary of the film - 4 


an increment in the film thickness At, = 0. 0375 R with 1 the radius of the sphere 
a being 1.1625 R. This leads to revision of the pressure diagram in Fig. 18(c), 
or rather of its base line, which should now be located at the level 0. 32(Ap) — .. 
from the pressure level 4, as determined by scaling of the distance KL mio | 
Fig. 18(a) and fitting it into the diagram, Fig. 18(c). 


new force of repulsion is found: 


“The exp expression for Pl is smaller than the earlier expression 
antes to a drier state of cement in which the one were located closer to- 


3 


the two relative positions of the grains considered here, he 


centres of the two adjacent grains hay have ‘moved er through the distance val q 
t. ty 0. 025 Re Ol R = 0.015 

At the same time the film thicknesses on both grains increased by the amount 

0.0375R each . The film decrease - - ratio is then 


————— = 5. This value of the ratio is comparable to the actual amounts 


Other in Shrinkage - - Swelling 
The of the shrinkage theory in this paper has been limited to a 


state of static equilibrium. Actually, however, the conditions are never static, 
_ since the cement is always either losing or gaining r moisture. The water film 
on the drying side of cement gr grain is thinner than on ‘the other side and in at- 
tempt to gain the equilibrium the water creeps or flows slowly along the 
a surface from the moist side to the dry side. The outer parts of the drying ff 
- specimen ; are always’ in a more advanced stage of drying and shrinkage than 
a _ the inner parts, and so the values of the weight and length of the specimen at 
_ all but the extreme moist and the extreme dry stages represent the mean 
; values ofa wide range of conditions. _ This is another reason why | the parr 
ratio earlier should be re- 


tad 


- thinks, hydration is never completed, and its progress is accompanied by a 
linear and volumetric changes. — No discussion of this important subject will 


this phenomenon is purely ch chemical, or perhaps, partly physical, as the author 
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- cause it ovdia strongly : at the surface while the centre ‘still ‘remains dry a 


_ therefore does not swell. The tension stresses thus set up in the interior — 
region may be greater than the ultimate strength and so the material cracks. a 
On the other hand, leaving a dry specimen in the open air causes a slower and— 
_more uniform absorption of moisture and swelling. A later immersion in ~~ : 
water completes the saturation. two- stepped procedure of this type is 
less severe as far as the integrity of the specimen is concerned than the anal a 
‘The most damaging operation from the viewpoint of cracking is the sim a: 
-mersion in water after the specimen is saturated with an inert liquid like a = 
= or oil. This is probably due to a particularly sharp division between oe 
the swollen and unswollen parts of the specimen, since the adsorption of rd 
- moisture takes place through the mechanism of the water film advancing slow- 
ly along the surface of cement grains. On the other hand, when a dry speci- _ 
men, devoid of inert liquid, is placed in water the transition between the © a g 
swollen and the unswollen regions is probably less abrupt. ft is believed that = : 
cracking of prismatic cement specimens of the type used in these studies can-_ 
not be completely avoided, although it can be reduced by proper precautions. b., 7 


‘Mortar specimens do not seem to crack on drying and possibly not on wetting. 


1 - Cracking of course is very damaging to the cement strength particularly _ 
ie the tension and the bending strengths. ‘Ast to its effect on the volume - — ¥ 
age graphs. it appears to be small. There are usually no more than one or vor 


two cracks on any face of the prism, and the crack on one face is almost never 

_ matched by a crack on another face at the same > section. Complete drying - 
usually ¢ closes the crack and makes it inconspicuous. The thickest crack ob- 
_ served was about 0.0007" but usually the cracks are smaller—0.0002" or ee 

0. -0003". . Compared ti to the total change in — of a 4" Specimen in the ae 


=A ing c curves of similar : specimens 5 subjected to ‘similar conditions. “One > would Fe — 
s think that appreciable variations in the shape of the graphs would be st 


_ The shrinkage studies were in to the earlier 
studies of creep in cement, and io have a direct bearing on creep. i 
_ As explained earlier, the moist cement grains are surrounded by films of ad- = 
_ sorbed water. The films are depressed at the point to point contacts of grains 
_in response to their intergranular attractions, delicately balanced against the 
_ repulsions produced by the film pressure. With water films in between, the 
— grains do not bear on each other directly. The film water, gy g under — 
‘| attraction of cement, is not a liquid but a plastic solid. some 
_ With this picture in mind it is easy to see what happens when a moist ae 
cement is to compression. On load, cement grains 


os 
ds to | | 2nd swelling. In this manner the author lost nearly all specimens made of ; * ae 
sphere| Cement finer than #325 mesh. _ Cracking is caused not by shrinkage and swell- a nei - 
-18(c), 
— | 
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r film _ the effect of cracking on the shape of the volume - shrinkage graphs seems 
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and water films de deform as solid elastic meee |The films directly between 
_the grains are so thin in relation to grains that they have practically no effect | 
on the elastic deformation. This explains the equality of the ‘moduli o of | elas- 

_ ticity and Poisson’s ratios for moist and dry cements during the rise or all : 
of load. (1) compression load increases the force tending to bring the 
grains closer together in the direction of the load. _ Before the application ae 
the load this force is A, - the grain to grain attraction. | After: the application 
of the load this force is A + G, the latter term being the effect of the applied = 

cae load. The force of repulsion of the grains, caused by the the film pressure, is 
still A, and consequently it is too small to balance the op opposing force, A + C. 

‘ef In order to re-establish the equilibrium the water in the film flows slowly as _ 

a plastic solid bring the grains closer together, while raising | the attraction 
: force to a greater value Ay and the repulsion to a still greater value Ay + +C. 
‘This slow adjustment of cement films isthe creep. 


" While the creep phenomenon thus agrees well in general with h the shrinkage 
§ of still require further study me their explanation. 


- All preceeding discussion is based on observations of shrinkage and swell- 

7 .: of cement prisms 3/4" x 3/4" x 4" subjected to various conditions of dry- 

ing and wetting, and described by means of graphs whose abscissae represent 

the internal volumes emptied or filled by the liquid, and the ordinates—the a 
Sa in length of prisms. The shapes of these graphs, and 

_ especially the behaviour of prisms after placing them first in inert liquid and 

later in water, ‘suggest th the explanation for the volumetric changes ir in the form 


- of the active water theory, according to which the water films are responsible 


for shrinkage, swelling and creep in cement. A detailed Giecussion of me- 


f chanical action of film water in the light of the theory is given. _ The theory is 
’ vr _ confirmed by the experimental evidence related to the behaviour of the drops | 
of water and kerosene on the cement surface and the rise of these liquids in 


‘The volume - shrinkage raphe contain a wealth of information still 
unutilized and even undigested. represent the outward manifestations 
of the physical and mechanical properties of cement and of the inner behaviour 
of this material under a variety of conditions. A change in the fineness and 
q kind of cement, water = cement ratio, conditions of hydration and drying, a age e, 
7 7 saturation with various inert liquids and other factors are immediately re- _ 


which to deduce the basic nature of the material and the laws governing it. | 


; - flected in the shape of the curves, giving information to the experimenter trom, 
+ 
— 


a _ Some peculiarities of the volume - shrinkage graphs still remain mysteries 


requiring further study, ott others been discussion 


The experimental work described here has been conducted in the Materials — é 


4 Testing Laboratory, Department of Civil Engineering, the University of a 
_ British Columbia. The work has bees ¢ given pumeaee financial assistance by 
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} the Department of Building Research, “The National Research Council of | 

— Canada. The author wishes to express his deep appreciation for this a 
sistance to the President of the National Research Council, Dr. ‘EL oe = 
Steacie, to the Director r of Building Research, Dr. R. F. ‘Legget, and to the 

4 Assistant Director of Building Research, Dr. N. B. Hutcheon. re 

: _ The e author also wishes to thank Mr. J. . F. Muir, Head, Department of Civil | 
Engineering, The University of British Columbia for his cooperation in al 
ing —— to the author the laboratory and id shop facilities of the easicesmmaall 


C. Powers and T. L. Bulletin 22, Research Laboratories 


the Portland Cement Association, 


4. Chemistry of Portland Cement by R. H. Bogue, 2nd edition, Reinhold 


Concrete Information #ST4, Portland Cement Association. — 


6 by A. H Hrennikott: Hydraulic Pressure in Concrete” by 


* Fundamentals of Hydro- yr Aeromechanics by Prandtl & uaa 


9. Elementary Mechanics of Fluids by Hunter Rouse, & 1946, >. 
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EC BCCENTRICALLY HINGED STEEL COLUMNS* 


‘by Theodore R. Higgins 


_ THEODORE 1} F. ASCE. —1 The test data presented by the 
a authors add further evidence that remarkably accurate predictions of column © 
behavior are feasible when—and only when —all of the factors influencing this 
behavior have been accurately estimated, and d provided the member is braced 
so as to eliminate premature lateral- torsional buckling. 3 Of greater interest, an 
perhaps, , than the new data are the authors conclusions, drawn largely 
‘Unless carefully read, the statement that the simple, widely used 
‘tion formula of Type I is unconservative can lead to misunderstanding. Ob- Bp -' 
viously, the type of formula can yield unconservative or grossly over- ae = 
conservative results, depending upon the values chosen for P’ and M'. The _ 
-AISC Specification for the Design, Fabrication and Erection of Structural _ 
Steel for Buildings is one of the several widely used design codes containing © = | 
the type; so far | as the writer is aware, , the : authors’ equation (I) is not cur- 
_ In Table 6 the data presented in Table 3 have been revised in terms of ata 
_AISC interaction formula, and also in terms of the same formula modified - 
_ injecting the factor (1 - 
3 Table: be compared with those given in Table 
‘Range’ of “Arithmetic Standard 
 Obse rvations 


"Interaction Type I 926 to 1. 0707 6.68% 
Typell 1. 1.117 te to 1. 384 0632 5.12% 
To eliminate | yield point as a factor in values F, corres- 
: - ponding to the slenderness ratios used in the tests, the following expressions 
were substituted the AISC ‘formula, is based upon a 


a. Proc. Paper Paper 1792, ‘October, 1958, by Richard E. E. Gordon | Fisher 


ss Director of irsis and F Research, American Ins Inst. of Steel Construction, New nen 


@ 
? 
— 
— 
; 
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Jys = 


INTERACTION FORMULAS - AISC WORKING STRESS RECOMMENDATIONS ~ 


690.250.9386 0.166 «1.102981 
3 69 0.75, 352 «1. 

3 «108 (0.75 0.715 0. 20 0 
= 1.50 0.506 0.600 1.106 

(0.988 0.2051 0. 
= 0.945 0.171 
5 0.711 O 
110.5 0.75 0.688 926 0. 
«1.50 0.468 0.510 
0.605 0.339 


+ 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 
x 


0.25 0.955 0. 162 1 
= 956 105 1.061 O. 
ae 0.595 0. 293 0. 988 
given in (1) in Table 6 were obtained by dividing the test 
i value for Pylt listed in Table 2 by 1.65 A. F,. As used in Table 3, M‘isthe | 
howe hinge moment, To convert the values for M/M' given in Table 3 to z, 
those comparable with the fp/Fp term in the AISC interaction formula it is 
only necessary to > multiply them by the the ratio of plastic bending mod modulus to 


q tie: It is seen that, taken in context with the rest of the code of which it is cur- # 5 
a part, the simple formula predicts the test results r re- 


— 
— 
| | ~ 
— 
| 
| | 
a 
| 
| 
¥ 
a 
| 
Bijlaard method, and the authors’ equation (II). On the other hand, when 
Type formulaerrors | 
the AISC working stress provisions are placed in the Typ formula errors 
lo ) ranging from 12 to 38 percent on the wasteful side are indicated by the same | 


bars, ars, unannealed 
I- Sections, unannealed 


Square bars, 
e - Sections, s, annealed alt 


~ 


is 


Bat: 


, in Fig. 9, 
: . results presented in Fig. 7 for the 22 end-restrained columns whose plastic 
bending modulus was equal to 1.5 times their section modulus. Based upon 
the data presented, there would appear to be little to recommend the injection — ee. 


the 1/(1 - P/PE) factor into the AISC Specification, unless at the same 
ti 


| 
ric Columns a 
= 
4] 
nearly sixty 


‘specitication(3) prior to the earliest adoption | of the AISC Specification. 
applied to compression members subject to combined axial and direct stress, 2 | Wi 
this formula was given as ; 


fers 

‘where S, the working stress in compression, applicable to both the axial and j y 
the bending term, was not permitted to exceed 16,000 - 70L/r. 
‘That this earlier design provision fell into disuse can be explained by the / a ~ 


unwarranted complexity of the de denominator of the bending term, The pre- 


“cision : implicit in the factor (1 - P/Pg) was incongruous with that of the work-— - 
ing stress provisions and with the overall accuracy of the analysis by which 
values for M and P are generally determined. Hence, resort to a simpler 
expression, | protected by a generous factor of safety, proved more popular 
__ The development of functions capable of predicting the strength of lateral- 
% ly and torsionally braced beam-columns of known end- restraint, subject to a 
a direct load alone and to bending alone, with greater accuracy than heretofore, - 
has again focused attention on the factor (1 - P/PE). However, specification 
provisions based upon these refinements will only place an added burden wll _ th 
; ce the designer, without commensurate improvement in the balance between _ P 
a safety and economy, iene complemented 4 an equally precise evaluation | of g 


and erected in compliance with them the past of 
century, have proven remarkably safe. tl 
- It is of interest to note that only through recognition « of the plastic rather ._ 


havior, were the authors able to achieve close correlation between computed 
_ and observed column strength. _ By the extension of this observation the use 
ii a of an allowable bending stress, higher than that permitted for axial tension, — 
also be providing the resulting deformations are 


ny 


— 

3 

— 

would certainly warrant 2 re-study of necessary eafety factor Des. 

| 

| 

2. The Theory and Practice of Modern Fra th Edition 
(1901) by J. B. Johnson, C. W. Bryan a bhn Wiley & <a 

The Design of St 
Ketchum - McGr 
— 


a eccentric cobaniha’ of a practically important type not hitherto investigated, 
_and to compare these data with various theoretical methods of strength pre- 


- diction. _ They have refrained from discussing the possible applications of a 
these findings to any particular conten codes, since this is not the purpose of 


a research report. im- 

; "plications « of these tests in n regard to the AISC Specification, are most wel i- 
- come as an indication of the practical relevance of this investigation. = == 
.. ‘Mr. Higgins iniaced to demonstrate that the AISC interaction formula of | 


bse II formula is ; overly conservative when applied to these tests. This 
demonstration | is achieved, however, , by utilizing a safety factor of 1 -65 in a 


for this formula is easily demonstrated to be not 1.65 but 2.0 (within a range ’ 
< of a few percent depending on L/r). 4 It is now established, and so stated a: Ls 
_ the AISC, that for beams the safety factor in that specification is 1.85 (see 
§ Plastic Design in Steel, Am. Inst. of Steel Construction, 1959, p. 2). The 

_ Specification conforms with general practice in | maintaining | a slightly ghar 


4 Higgins, gives unconservative results when applied to the present tests. “In 
fact this is so for 20 of the 24 eccentric ‘column tests, with deviations on the - 
It is evident that for an interaction formula to be consistent, the terms P : 
and M' cannot be arbitrarily “chosen” "as implied by the discussor, but 
= to the best of available knowledge, the calculated strength of _ — 
‘particular member when loaded as a concentric column on the one hand (P’), | 7» a 
and as a flexural member on the other (M' ). For the hinged specimens at s 
bent it is well established that M' must be the plastic moment of the section, © 
; and P' either the Euler value for long columns or the Engesser- “Shanley value 


(see. Column Research : “The Basic Column Formula,” Techn. Memo. 
No. % & 1952). These values have been so used by the writers, and Mr. ‘Higgins’ 


However, by using P' = 1.65-A- F, as the denominator of the first 
Mr. Higgins maintains that the actual strength of a concentric, hinged, pe 


| 


ft _ Asst. . Prof., Dept. of Structural Engineering., Cornell Univ., Ithaca, N.Y. ‘ 
Associate Prof., Dept. _of Structural Eng., Cornell Univ., Ithaca, 


| RICHARD E. MASON,! GORDON P. FISHER,” F. ASCE AND GEORGE 
q 

“ 
— 
ramers of this specification would be the first to object ifit were implied 
— the reverse is true. As will be shown below, if the correct safety factor eS ei = 
— 0 is used in connection with the AISC column formula, then the AISC in- [——— 
i 
‘2 
— 
3 


— = | is 
— 

| 


is: entiatactorily by multiplying AISC by 1.65. ‘It is s easily 
demonstrated that this is not so. A very short concentric column see il 
| a fails by simple yielding. Since the specified minimum yield em 
for ASTM A-7 steel is 33 ksi and the AISC stress for 17 


ze so that the latter must be rmultiplied by 20 20.56/10.02 = 2.05 (rather than 
A 65) to result in the correct strength, This, is an abbreviated way, demon- — 
_ strates that, with minor deviations, ‘the load factor implied in the AISC _ 


‘The correctness of this fact is demonstrated further by evaluating the a : 


_ concentric columns tests reported in Table 4 of the paper. For these pc 

_ types of specimens the actual yield points were within 5% of each other, om 
that it is justified to take the average yield point, 43.5 ksi, for all six tests. 7 
For this yield point the AI AISC f formula, modified exactly as as was done > by Mr. tdi 


This « expression is bracketted, as it must his formulas for the two 
€ pertinent steels with 42.5 and 44.5 ksi, respectively. Figure 10 shows the i 
k results of the six concentric tests of Table 4. Curve (1) of Fig. 10 shows the e y 
e above AISC formula as used by Mr. Higgins, i.e. », multiplied by 1.65. é It is fees} 
seen that this strength prediction is completely at variance with the tests and 
underestimates the real column capacities by a wide margin. . Therefore, 


_ cannot be u used for P' i in the interaction formula, On the other hand, curve (2) 


2.0. Tt is seen thi that this curve ; agrees” reasonably well with the concentric t test a! 


results even though it implies considerably larger residual stresses than \aeglae 
were present in the specimens, Also shown is curve (3) for o .. and d_,as _ 
used by the writers for determining Pp. It is seen to be in satisfactory agree-- 


|= with these concentric tests although it is evident that this curve, — 


wis, is based on larger residual stresses than were actually present. * 4 
= noted in the paper and accounts for the deviation of ¢ cES on the low side a 

the ‘range of moderate L/r- ratios. 

i _ Having thus established, by theory and by comparison with tests, that the 


load factor applicable to the AISC column formula is 2.0 rather than 1.65, the | 
writers have recalculated Mr. Higgins’ Table 6 on this basis. The results are aa 


formula of Type I overentinelée the strength of all but four of the 24 test 
columns, by margins up to 20%, | On the other hand, the Type 0 formula a 
conservative in all cases, in many instances by a considerable margin which = 
is due to the excessively low proportional li: limit implied in the saesnimmedl 

formula (when applied to these tests). 


Evaluating Table 7 in the same manner as the discussor did Table 6, 


— 

.4 

5 

. 


Interaction 0,806 to 1. 020 0.9 0. 0655 7. .06%, 


Interaction 1,017 to 1. ok 100 
‘Type I 


- corrections for this set of tests, since their result is evidently of the same 


nature as just demonstrated for the hinged columns. 


-_ __ While the writers believe to have shown that Mr. Higgins’ re-evaluation of a Pi 
their tests is distorted by an oversight in regard to the true safety factorof 


the relevant AISC formula, they appreciate Mr. Higgins’ attention to these _ ae 
tests and the thoughtfulness expended on their 


at 
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Ty raphical errors i in original pape er: 
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8 Coefficient of term of expression for should be be 5, 000° 


nstead o 


9 first line, inequality should read 


Pa e 10 line 14 from bottom. _ Expression in denominator of second term 
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Fig, 9, concerning the tests o 
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DISCUSSION 


iy MR. » HIGGINS" TABLE 6 CORRECTED FOR AISC COLUMN FACTOR OF OF SAFETY = 2.0 * 


0.25 0.772 «0.186 


0.166 


a 76 a3 0.889 06435, 


75 0.590 0.260 0. 0. 184 


0. 0.815 0.205 “2.020 0.2h2 
0. 25 0. 951 0. 1.082 


0.75 0. 568 0.238 0.806 806 1-105 

36 a 1. 250 | 0. 0. 

0.649 1.035 
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Discussion by Enzo Oscar Macagno. 


“te 


_ For a long time and in many publications the rule mentioned by the authors: ‘4 
in Step 2 has been given. The writer wishes to point out that although the 
rule works in most of the cases, it is not ¢ correct and fails sometimes. For — 


some wise remarks about dimensional analysis. 


of variables is 3 equal to the number | of dimensions. In fact, we know the ar an- . 
; swer is given byc/VE/P as the dimensionless group. The correct rule 
. based on r being the rank of the saa eeanes of the number of dimen- 


sions involved, 


Proc, January, 1959, by Charles Cc. and Vaughn E, 
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4 = DIVISION ACTIVITIES 


ENGINEERING ME MECI CHANICS DIVISION 


of the. of Civil | Engineers 


AT WASHINGTON, C. CONVENTION 


‘iia Professor E . Popov has announced the program for three half day mibiens 


at the October 19- 23, 1959, ASCE Annual Convention sponsored ey by = 


_ the Engineering Mechanics Division as well as six half day sessions co- ae on 


Monday Morning—October 19—Column Research 


a 
ei (1) “Column Research Council”, B. G. Johnston, Univ. of Michigan 
4 (2) “Basic Column Strength”, L. S. Beedle Univ. 
(3) “Effective of Framed Kavanagh, New York City. 
Monday Afternoon—October 19~column of 


capers with the Structural Division. The for these sessions are 


“Buckling Problems in in Plate Girders*, Basler and B. “Thurlimann, 


Note: No. 1959-27 is part of the copyrighted Journal of the Mechanics 


Proceedings of of the American a Society of of Civil Engineers, Vol. 85, EM: 


. __—* 1959 by the / the American Society of Civil Engineers. 


a 
— 
| 
‘ 
“Lateral Buckling of Beams”, J. W. Clark and H. N. Hill, Aluminum Co. 
Ake (5) “Postbuckling Strength and Effective Width of PlatesinEndCom- 
4 
— 


Tuesday Afternoon—October 2 20—Column Research 


(10) “The Effect of Floor System Participation on Top Chord Stresses of a . 
Single | Span Pony Truss R. M. Barnoff, 


_ Wednesday - Fluid of Stratified 


; 
(1) Int Introductory Remarks, D. R. F. Harleman, M.I.T. es eo 


(2) “Stratified Flow Into a Line Sink for a . Linear Density Gradient”, 
= Debler, Univ. of Michigan. 
(3), “Irrotational Motion of Two Fluid Strata Towards a. Line sink’, 
a vd G. Huber, McMaster  Univ., Ontario. $ se 
“Internal Waves in a Stratified Fluid—A General Review”, G. H. 
Morning—October 22— —Lightweight Construction in Metals 


(1) “Welded Aluminum Bridges for Military Traffic”, J. Tamanini and 


“Tests of a a Composite Aluminum and Concrete Highway Bridge” 


a (3) “Design of Welded Aluminum Structures”, H. N. “HU, J. Ww. Clark, , and 
Brungraber, Aluminum Co. of America, 
4) “Strength of Welded Aluminum Columns’, R. J. Brungraber and #. 
J. W. Clark, Aluminum Co. of America. 
hursday Afternoon—October 22— Construction in 
6) “Light Gage, Cold- Formed Steel Construction’, G. Cornell 


a (6) “Stability of Thin-Gage Metal Structures”, Uni 
“Postbuckling Behavior of Rectangular Plates” , M. ‘Stein, 


: as Large Deflection Vibration of Plates”, G. Herrmann, Columbia Univ. - 


(Thursday sessions co-sponsored with Structural Division.) _ 
- Friday Morning—October 2 23—Structural Dynamics 
— R. J. Hansen, » presiding ites 
; — “Dynamic Elastoplastic Response of Rigid Fr Frames” »] F. L DiMaggio 
“Dynamic Effect of a Moving | Load on on a Rigid Frame”, 
“Effect of End Fixity on the Vibration at Rods” »D. Nuclear 


| 
| 

= 

gt 
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Friday Afternoon—October 23—Analysis of Highly Redundant Structures by 


| ‘Experimental and Analogue Techniques. J. 8. Archer, 


(1) “Electronic Models for Structural Vibrations with Elastic-Plastic 

Couplings”, H. M. Paynter, American Center for Analogue Computing. 

(2) “Aircraft Structural Analysis with the Electric Analogue Computer”, 
w. Brignac and R. G. Schwendler, Convair, Fort Worth. ts. 

a (3) “Arch Dam Analysis with the Electric Analogue Computer”, _R. H. Mi 

=) “Elastic Model Design of the B-58 Airplane Structure”, J. 'W. Wells and a ; 


= 


B. Convair , Fort Worth | 


ae ‘Four papers at this conference (Stanford, California, September 9-11) are 
‘gponsored by the Engineering Mechanics Division of ASCE. The complete 
4 program will be arranged from these papers (see final selection below) | <a * ; 


(1) “Limit Analysis of Simply Supported Circular Shell noaanll 
Fialkow, U.S. Army,New York, 
(2) “Large Deflection of Elasto- Plastic Plates U ‘Under Uniform Pressure”, 
___ Thein Wah, Southwest Research Institute. 
sc (3) “On the Asymmetrical Bending of Conical Shells”, B. Wilson, ey of a4 


California, Be Berkeley, and M. ant: M. L. Williams, California Institute 
INTERNATIONAL SYMPOSIUM | ON FLUID IN THE 
Professor R. E. Fadum calls the divisions attention to the above s iladie. 
um being held at Cornell University, July 8-15, 1959. Its aim is to bring to- 
ome international experts in the fields of Fluid Mechanics and lonospheric | 
Physics with the intent of gaining a more complete and satisfactory under- — 
; - standing of the dynamical processes occurring in the ionosphere. Emphasis = 
__ will be placed upon the phenomenon of turbulence and its effects. ae ee 
= The meeting is being sponsored by the International Scientific Radio Union — se 
BS. _ together with the International Union of Theoretical and Applied Mechanics, yee : 
the International Union of Geodesy and Geophysics, and the International 
~ Astronomical Union; and with support from UNESCO and the U. S. ‘National yan" 
ie a Be The major part of the Symposium will be the sessions on Monday, Tuesda > 
e and Wednesday, July 13-15. July 8-10 will be devoted to semi-formal presen- 
_ tations for the purpose of exchanging existing knowledge between the fields a 
_ Fluid Mechanics and Physics. 
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COLUMN PAPERS AVAILABLE 


Column Research Council of the Engineering Foundation has released 


ie following list of publications currently available from the Secretary of acme * 
> - Research Council. Requests may be addressed to: Secretary, Column Re- 
= Search Council , 319 Engineering University of Michigan, Ann Arbor, 


- “The Basic Column Formula”—Technical Memorandum No. 1 


“Design Charts for Elastically Restrained | Eccentrically Loaded Columns”, 


Determination of the Buckling Load for Columns of Stiffness”, 


7 sat of Initial Eccentricities on Column Performance and Capacity”, 


“German Buckling Specifications” —DIN 4114, Vol. 1—German Version 
“German 4114, ‘Vol. Version 


“Japa anese »se Handbook | on Elastic Stabilit (in Jar inane)” 
“Lateral Buckling of Eccentrically Loaded I- and H-Section Columns” 


“Notes on 1 Compression Testing of Me Metals” '—Technical Memorandum No. 


The Philosophy of Column 2”—Column Council 


“The Stability of Bridge Chords Without Lateral Bracing”, R. M. Barnoff 


Study of Columns with Perforated Cover Plates”, White an 


“A Survey of Progress 1944-1951”, B. G. Johnston and the Project rh 


“Tests Analysis of Eccentrically Loaded Columns*, R. E. Mason, 
BARLY TRANSACTIONS VOLUMES OBTAINABLE 
ae feasibility of reproducing the first ten volumes of ASCE Transactions ~ Sho 
2 —_ -1881) has been studied. It has been decided that these historic volumes __ 
g ‘fe could be reproduced at a cost that would permit a top price of $150 for the ae 
in obtaining such a set, the p1 project will be undertaken. - it the endeavor is 
Es ——— other rare volumes of Transactions will be reprinted. = 
Engineers interested in obtaining the ten-volume set should write tothe 
Executive Secretary of ASCE, 33 West 39th Street, New York 18, N. Y. 
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